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ABSTRACT 


t 


Atmospheric  turbulence  data  in  many  forms,  from  many  sources,  have 
been  combined  in  both  spectral  and  exceedance  probability  models  of  the 
vertical  component  of  atmospheric  turbulence.  These  models  will  be  useful 
in  aircraft  structural  design  procedures  as  well  as  in  other  analytical 
procedures  involving  aircraft  operations. 

For  the  spectral  probability  models,  a  simple  analytical  form  for  the 
spectrum  is  chosen  which  yields  a  good  fit  to  the  available  measured 
spectra.  The  parameters  of  the  spectrum  are  related  to  physical  and 
meteorological  factors  in  tabular  and  graphical  form,  and  probability 
distributions  of  the  pertinent  meteorological  parameters  have  been  obtained 
from  climatological  data.  The  resulting  probability  distributions  of  the 
spectral  parameters  have  been  computed. 


A  normalized  exceedance  model  is  developed  that  represents  the 
overall  cumulative  probability  distribution  of  the  vertical  turbulence  velocity 
for  the  special  probability  density  function  of  the  rms  vertical  gust  velocity 
recommended  in  NACA  TN  4332.  The  basic  model  involves  one  parameter, 
namely  altitude,  but  may  be  modified  by  using  correction  factors  that  have 
also  developed.  The  correction  factors  account  for  varying  terrain  types, 
season,  and  time  of  day. 

These  models  of  the  vertical  component  of  atmospheric  turbulence 
represent  what  should  be  an  improvement  over  previous  models  because  of 
the  significant  amount  of  data  which  has  become  in  recent  years.  Never¬ 
theless,  the  lack  of  sufficiently  reliable  data  remains  a  substantial  problem. 
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SECTION  1 
INTRODUCTION 


One  of  the  primary  problem  areas  in  the  design  of  airplane  structures 
is  the  atmospheric  turbulence  environment.  Turbulence  may  be  the  critical 
structural  design  problem  through  either  the  small  repetitive  gusts  or  the 
large,  infrequent,  but  sometimes  disastrous  gusts.  Moreover,  for  some 
airplanes,  particularly  transport  airplanes,  both  the  large  and  the  small 
gusts  can  lead  to  critical  design  problems.  The  large  gusts  determine  the 
basic  structural  static  strength,  and  the  small  gusts  determine  the  fatigue 
life  of  the  structure.  To  design  properly  the  airplane  structure  to  withstand 
the  rigors  of  the  atmospheric  turbulence  environment  within  acceptable 
reliability  limits,  a  model  of  the  turbulence  environment  must  be  available 
in  a  form  amenable  to  the  modern  techniques  of  structural  load  calculation. 

The  statistical  methods  now  in  use  appear  to  be  the  most  logical  approach 
in  view  of  the  statistical  nature  of  the  turbulence.  What  is  needed,  then, 
is  an  adequate  and  reliable  model  of  atmospheric  turbulence  that  is  amenable 
to  the  statistical  methods  for  load  calculation. 

The  work  covered  in  this  report  represents  an  effort  to  update  and 
improve  the  models  of  the  vertical  turbulence  velocity  now  in  use  by 
including  new  and  more  recent  data  as  well  as  by  trying  new  analytical  methods. 

A  significant  quantity  of  new  spectral  and  VGH-type  data  has  become 
available  since  publication  of  the  NACA  (NASA)  TN  4332  in  1958  (Reference  27) 
and  MIL-A-8866( ASG)  in  I960  (Reference  53),  the  latter  being  based  strongly 
cn  the  results  of  the  former.  Relatively  large  quantities  of  spectral  and 
VGH  data  for  altitudes  below  1000  feet  above  the  terrain  have  become  avail¬ 
able  recently. 

Furthermore,  a  large  percentage  of  the  spectral  and  some  of  the  VGH 
data  were  presenter  with  enough  accompanying  meteorological  and  physical 
data  to  permit  a  more  sophisticated  model  of  turbulence  in  the  atmospheric 
boundary  layer,  Through  recent  sampling  programs  employing  the  U-2 
airplane,  a  largu  quantity  of  VGH  peak  count  data  has  become  available  for 
altitudes  approaching  75,000  feet.  These  data  are  especially  concentrated 
in  the  altitude  band  from  50,  000  to  70,  000  feet  where  other  aircraft  have 
not  been  able  to  operate  satisfactorily  for  sampling  purposes,  On  the  other 
hand,  except  for  spectra  measured  during  thunderstorm  and  cumulus  cloud 
penetrations,  there  are  no  spectral  data  available  for  altitudes  above  some 
very  few  thousand  feet. 
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The  turbulence  modeling  problem  has  been  approached  from  two 
essentially  different  directions: 

1,  An  analytical  form  for  the  spectrum  is  chosen  so  that  the 
relations  between  the  spectral  parameters,  scale  and  mean- 
square  intensity,  and  the  meteorological  parameters,  primarily 
stability,  mean  wind  speed  and  height  above  the  terrain,  can  be 
identified  and  determined.  Then  probability  density  functions 

for  the  meteorological  parameters  are  determined  from  climatolog¬ 
ical  data. 

2.  Based  on  the  assumption  of  a  spectral  function  including  the 
scale,  but  not  the  mean  square  intensity,  measured  peak  count 
and  exceedance  statistics  from  normal  acceleration  data 
recorded  during  routine  commercial  and  military  flight  operations 
are  converted  into  both  the  probability  density  anu  the  cumulative 
probability'  of  the  root -mean-square  intensity  of  turbulence. 

Both  approaches  have  been  designed  and  carried  out  to  provide  models  of 
the  actual  atmospheric  turbulence  environment,  models  that  are  presumed 
to  be  independent  of  any  given  airplane,  airplane  operations,  geographical 
locations,  and  other  sources  of  unwanted  bias.  The  two  approaches  were 
employed  simply  because  of  the  basic  theories  relating  measured  data  on 
the  one  hand  and  the  parameters  with  which  the  measured  data  vary  on  the 
other.  The  spectral  data  are  most  easily  related  to  meteorological  param¬ 
eters.  The  peak  count  or  exceedance  data  are  most  easily  related  to  those 
parameters  associated  with  airplane  operations,  especially  altitude.  In 
the  range  of  altitudes  near  the  ground  where  both  spectral  and  peak  count 
or  exceedance  data  abound,  it  should  be  possible  to  obtain  good  agreement 
between  the  results  of  both  approaches,  at  least  within  a  restricted  set  of 
conditions  that  would  in  theory  permit  agreement. 
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SECTION  2 
SUMMARY 


As  discussed  in  the  Introduction,  the  atmospheric  turbulence  model 
is  developed  in  this  report  along  two  essentially  different  approaches,  one 
involving  the  turbulence  spectra  directly  and  the  other  involving  peak  count 
and  exceedance  statistics  in  which  the  spectra  play  an  important  but  indirect 
role.  The  different  approaches  were  employed  because  of  the  nature  of  the 
basic  data  and  because  of  the  analysis  procedures  that  were  employed  to 
develop  the  turbulence  model,  The  model  was  developed  only  for  the  vertical 
component  of  the  turbulence  velocity. 

Section  3  deals  with  the  turbulence  power  spectra  and  their  relation¬ 
ship  to  the  physical  and  meteorological  parameters: 

1,  surface  roughness  and  terrain  type, 

2,  height  above  the  surface, 

3.  wind  conditions  (both  velocity  and  shear),  and 

4.  atmospheric  stability. 

An  analytical  form  for  the  spectrum,  that  proposed  by  Lappe  and  Davidson 
(Reference  1)  and  later  used  by  Lappe  in  his  analysis  of  B-66  low-level 
gust  data,  was  chosen,  It  appears  to  fit  all  the  available  measured  spectra 
within  limits,  yet  is  simple  enough  to  facilitate  straightforward  analysis. 

This  spectral  form  fits  not  only  the  B-66  low-level  data,  but  also  the 
available  tower  spectra,  low-level  spectra  measured  by  aircraft  in  England 
and  Africa,  and  thunderstorm  penetration  spectra. 

The  relationships  found  between  the  spectral  parameters  and  the 
physical  and  meteorological  factors  are  presented  in  tabular  and  graphical 
form  suitable  for  design  and  prediction  purposes.  The  distributions  of 
the  pertinent  meteorological  parameters  have  been  obtained  through  a 
climatological  study,  and  the  resulting  expected  distributions  of  the  spectral 
parameters  computed  using  the  proposed  interrelationships. 

Thunderstorm  spectra  are  analyzed  separately  in  Section  4  because, 
unfortunately,  there  are  little  if  any  supplementary  data  to  accompany  the 
measured  spectra.  It  was  found  difficult  to  correlate  measured  spectra 
with  any  physical  or  meteorological  parameters.  The  B-66  storm  penetration 
spectra  were,  however,  fitted  with  the  previously-mentioned  spectral  function, 
and  it  was  found  that  a  possiVv  useful  correlation  apparently  exists  between 
the  mean-square  intensity  (varicc;e)  and  the  scale  length.  The  scale  length 
increases  with  Increasing  metn^iquare  intensity. 
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Section  5  is  devoted  to  a  brief  summary  of  the  existing  knowledge  of 
clear  air  turbulence  and  related  physical  and  meteorological  parameters 
at  the  higher  altitudes,  for  which  at  present  there  are  little  data  of  any 
statistical  significance. 

Through  the  efforts  of  NASA  and  various  military  groups  in  the 
USA,  England  and  Canada,  a  considerable  quantity  of  VGH-type  statistical 
data  has  become  available  in  the  past  few  years.  Fortunately  these  more 
recent  data  have  been  mostly  well  documented  and  well  supported  by 
the  necessary  supplementary  information  needed  to  permit  accurate  analysis 
of  the  data  and  to  make  reliable  estimates  of  any  statistical  bias.  The  VGH 
data,  along  with  certain  spectral  data  that  existed  in  statistically  signifi¬ 
cant  quantities,  were  developed  into  an  exceedance  probability  model  of 
turbulence.  The  particular  graphical  form  of  the  model  is  the  form  recommended 
by  Houbolt,  Steiner  and  Pratt  (Reference  24),  that  is,  the  results  are 
presented  as  "generalized  prediction  curves". 


In  particular,  the  exceedance  model  consists  of  parameters  defining 
both  the  probability  of  encountering  turbulence  of  a  given  type  and  the  intensity 
of  the  given  type  of  turbulence.  These  parameters  are  presented  primarily 
as  a  function  of  height  above  the  terrain.  Correction  factors  are  given,  when 
possible,  to  account  for  the  variation  of  turbulence  with  physical  parameters 
other  than  height  such  as  terrain  roughness,  season  and  time  of  day.  The 
correction  factors  are  presented  generally  only  for  low  altitudes  above  the 
terrain  since  it  is  only  at  the  low  altitudes  where  sufficiently  reliable  data 
exist. 


The  exceedance  model  is  presented  for  all  altitudes  from  ground  level 
to  100,000  feet,  but  above  approximately  70,000  feet,  reasonable  extrapola¬ 
tions  were  made  because  no  directly  applicable  data  exist. 

Generally,  the  exceedance  model  presented  in  this  report  may  be 
compared  with  the  earlier  NASA  (NACA)  model  in  Reference  27  .  The 
comparison  will  show  that  for  small  vertical  gust  velocities  the  model  of  this 
report  is  more  severe  (i.  e.  ,  has  a  greater  number  of  exceedances),  whereas 
lor  the  larger  gust  velocities  the  present  model  is  considerably  less  severe. 

Although  the  spectral  and  exceedance  models  of  the  vertical  gust 
velocity  presented  in  this  report  are  considered  to  be  improvements  over  and 
extensions  of  previous  models,  the  resulting  models  clearly  show  the  effects 
of  the  lack  of  sufficient  data.  Obtaining  meaningful  and  useful  atmospheric 
turbulence  data  at  all  altitudes  of  interest  remains  the  area  in  which  the 
greatest  efforts  should  be  directed. 


4 


AFFDL-TR-65-122 


SECTION  3 

THE  LOW-LEVEL  TURBULENCE  MODEL 

i 


This  section  deals  with  turbulence  in  the  layer  of  atmosphere  extending 
from  the  surface  to  about  1C00  feet.  The  boundary  layer  occupies  approximately 
the  first  300  feet  of  this  depth  with  the  remainder  being  part  of  a  gradual  tran¬ 
sition  to  the  geostrophic  region  above.  It  is  in  this  1000-foot  layer  that  most 
of  the  turbulence  spectrum  measurements  have  been  obtained  and,  hence,  the 
greatest  potential  exists  for  fruitful  analysis. 

The  data  used  for  the  low-level  spectral  analysis  were  obtained  both 
from  fixed  tower  measurements  and  from  aircraft  measurements.  The  tower 
data  came  mainly  from  two  sources.  Round  Hill  (Reference  2),  and  Brookhaven 
(Reference  3).  The  aircraft  spectra  were  taken  from  the  Douglas  B-66  pro¬ 
gram  (Reference  4)  and  from  an  R.  A.  E.  low-level  study  (Reference  5).  If 
both  types  of  measurements  are  correct,  the  data  must  necessarily  fit 
together  in  a  fairly  orderly  fashion  to  produce  a  continuous  and  logical  model 
over  the  1000-foot  layer. 

The  results  of  the  analysis  do  indeed  show  that  the  tower  and  aircraft 
data  are  compatible.  Since  the  tower  data  have  much  more  detail  in  the 
accompanying  meteorological  measurements,  they  have  proven  valuable  in 
determining  the  meteorological  factors  affecting  the  higher  level  aircraft 
data  where  reliable  stability  measurements  do  not  exist. 


3.  1  THE  FORM  OF  THE  SPECTRUM 

Monin  (Reference  6)  has  shown,  using  similarity  theory,  that  the  spec¬ 
trum  near  the  ground  can  be  expressed  as 

f  §  (-f)  —  (F,  2/r)  ( 1 ) 

whe  re  <j(f,  z/r)  is  a  universal  function  and  f  $  (f)  is  the  logarithmic  spectrum, 
i.  c.  ,  the  product  of  frequency  and  the  spectral  density.  The  Monin-Obukhov 

stability  length  P  is  defined  by (u.*€  /°Tj/(yKH\  where  Z  ,  V  is 

the  surface  stress,  < ^  is  the  specific  heat  of  air  at  constant  pressure, 

K  -  0.  40  is  von  Karman's  constant,  and  p  is  the  vertical  turbulent  heat  flux. 
Equation  1  states  that  the  shape  and  magnitude  of  the  spectrum  is  determined 

completely  by  0"w  ,  f  ,  and  Z/p  .  The  similarity  theory  can  be  expected  to 

apply  only  at  relatively  low  levels  (below  100  meters,  perhaps)  and  over  rela¬ 
tively  smooth,  homogeneous  terrain.  However,  the  same  spectral  shape 
appears  valid  for  a  much  broader  range  of  heights  and  conditions.  Since  the 
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boundary  layer  case  must  eventually  fit  within  any  over-all  turbulence  model, 
it  is  not  unreasonable  to  begin  by  generalizing  the  dimensional  argument  pre¬ 
sented  in  Equation  1  to 


-  o£s(F.'*/n' 


(2) 


whe  re 


The  new  length  parameter  t.  replaces  2 


We  would  expect^ 


to  be  equal  to  2  rear  the  ground  over  smooth  terrain,  but  to  possibly  depart 
from  this  behavior  in  other  conditions.  There  is  some  inconsistency  in 
replacing  Z  by  2  in  F  ,  but  not  in  2  fp  ;  however,  at  this  point  we  carry  the 
2/p  dependence  only  as  a  reminder  that  the  spectrum  is  a  function  of  stability. 


From  a  general  review  of  the  literature,  we  are  led  to  consider  two 
possible  analytical  expressions  for  ( fj .  From  Panofsky  and  McCormick 

(Reference  7) 


f§Jf)  -- 


cr* 


(’+  4P)‘ 


(3) 


and  from  Lappe  (Reference  1) 

(itzn-Z-ff 


Both  of  these  models  have  proven  capable  of  fitting  observed  spectra.  The 
former  has  been  applied  mainly  to  data  from  fixed  towers;  the  latter  to 
aircraft  observations.  It  is  readily  seen  that  these  two  expressions  yield 
the  same  spectral  shape.  However,  neither  includes  the  effect  cf  stability, 
and  the  Lappe  equation  does  not  explicitly  include  any  influence  of  height, 
although  </  turns  out  to  be  a  function  of  height  when  Equation  4  is  applied 
to  observations. 


To  relate  the  variables  in  the  two  forms  so  that  a  single,  consistent 
model  can  be  applied  to  both  kinds  of  data,  and  to  bring  both  models  into 
agreement  with  similarity  theory,  we  generalize  the  Panofsky- Me  Co  rmick 
modei  by  writing 


cr  A  )  F' 


\jf(a  (*/r)F'J 

We  have  introduced  a  function  of  stability  A  (*/p) in  place  of  the  constant  4, 


(5) 


and  we  have  replaced  F  by 


We  now  have 


i(Wn) 


_ A(2/r)F' 

[>/  <4  (  Z/P)F  'J 


(6) 
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in  accord  with  similarity  theory.  It  follows  that,  for  the  Lappe  model  to  be 
consistent,  we  must  have 


A(i/r)t 


or 


(7) 


We  now  have  an  analytical  expression  for  the  spectral  density  which  we  may 
write  in  terms  of  the  longitudinal  wave  number  A  as 


(A)  = 


*  A  i 


(i+a/4)  1 


(8) 


where  A  is  a  function  of  stability  and  Jzziiox  low  levels. 

The  analytical  expression  of  Equation  8  has  been  checked  for  fit  against 
33  spectra  computed  from  tower  wind  measurements  and  over  100  spectra 
computed  from  low-level  gust  data  taken  by  aircraft.  The  method  used  to  fit 

the  curves  was  simply  to  solve  Equation  8  for  <T^  and  At  using  two  points  on 

the  measured  spectra.  The  two  points  are  chosen  near  the  upper  and  lower 
irequency  limits  of  each  speccruin,  but  far  enough  from  the  ends  to  insure 
that  the  data  are  free  of  errors  inherent  in  the  measurement  technique.  In 
every  case,  a  very  close  fit  can  be  obtained  by  suitable  choice  of  these  points. 


Transformation  of  spectral  parameters  is  made  in  the  following  way: 


cd=f tjndf  =  j §„Lf(*)idftA)  =  f 


1 

By  notation 


$  */?(*)/  zi  =  §<*) 


d  4 


A  X  4 

so  that  °us~  C  $  (f)Af  —  f  $  (A)  <AA 

J  y-  J  “S 

0  A, 


Now  for  the  particular  transformation  from  Equation  5  to  Equation  8,  we 
use  Taylor's  hypothesis,  f  -  AV ,  so 

$  (A)  -  v  i'f}-,  A=o, 
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If  we  differentiate  A  §JA)  with  respect  to  A  ,  we  get: 

z  [~  At  _  2  A  (Alf | 

Equating  this  expression  to  zero,  we  find  the  logarithmic  spectrum  has  a 
maximum  when 


(9) 


and  Jkm  is  the  wave  number  at  the  maximum  point.  Thus,  we  see  that 
corresponds  to  the  wavelength  of  maximum  energy  in  the  spectrum.  Using 
Equations  8  and  9,  it  is  possible  to  determine  AJL  and  cr ^  from  a  single  point 
(the  point  of  maximum  energy)  on  the  logarithmic  spectrum. 


3.2  THE  FIT  TO  THE  TOWER  SPECTRA 

The  results  of  the  Round  Hill  experiments  include  both  the  logarithmic 
spectra  and  the  spectral  density  curves.  At  these  lower  levels,  the  lower 
frequency  limit  of  the  measurements  is  well  below  the  peak  of  the  A  $ ^  (A) 

spectra  and  ti.e  peak,  in  most  cases,  is  well  defined.  These  data  afforded 
an  excellent  opportunity  for  checking  the  At  value  predicted  by  the  two-point 
solution  of  Equation  8.  The  two-point  solution  was  obtained  using  the  spectral 
density  curves  and  the  value  of  Ai  obtained  from  this  solution  was  compared 
with  the  value  determined  by  Equation  9  and  direct  observation  of  Am  on  the 
§J*)  spectra.  The  agreement  is  very  good,  indicating  that  the  shape  of  the 
spectrum  near  the  ground,  ever,  at  the  long  wavelengths,  is  closely  represented 
by  the  proposed  ana lytical  model.  Figure  1  shows  some  representative  com¬ 
parisons  of  measured  spectra  and  those  computed  by  a  two-point  fit  of  the 
associated  spectral  density  data  to  the  proposed  analytical  equation. 

3.  3  THE  FIT  TO  AIRCRAFT  SPECTRA 

In  general,  the  aircraft  data  do  not  cover  low  enough  frequencies  to  be 
able  to  measure  •  However,  there  has  been  one  set  of  aircraft  measure¬ 
ments  made  (Reference  8)  that  covers  a  range  of  wavelengths  from  1  0  ft  to 
60,000  ft.  A  composite  spectrum  from  this  work  is  shown  in  Figure  2.  The 
superimposed  dotted  curve  is  the  computed  spectrum  obtained  by  fitting 

Equation  8  to  the  points  A  --.  0005,  and  A  -  .05.  The  agreement  is  exception¬ 
ally  good  over  the  whole  wavelength  range,  including  the  low-frequency  end. 
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k  CYCLES/FT 


Figure  I.  Comparison  of  Measured  Tower  Spectra  w'th 
Analytical  Form 


.00001  .0001  .001  .01  .1  | 

INVERSE  WAVELENGTH,  CYCLES/FT 

Figure  2.  Comparison  of  Aircraft  Spectrum  with  Analytical  Form 


AFFDL-TR-65-  122 


Lappe  has  shown  that  the  chosen  analytical  eauation  fits  the  B-66  spectra 
over  the  measured  wavelength  range  and  the  work  done  at  CAL  with  these 
spectra  certainly  confirms  this.  A  very  good  fit  lias  also  been  obtained  to 
spectra  published  by  the  Royal  Aircraft  Establishment  in  England  (Refe rence  5). 

It  has  been  well  established,  during  work  on  this  contract,  that  the 
spectral  density  Equation  8  provides  a  good  fit  to  all  the  observed  spectra 
over  the  measured  range  of  wavelengths.  From  the  tower  spectra  and  from 
the  spectrum  measured  by  Crane  and  Chilton  (Reference  8),  it  appears  that 
this  spectral  form  provides  a  realistic  fit  to  the  actual  spectral  density  down 
to  frequencies  at  least  below  the  peak  energy  point. 


3.4  TURBULENCE  INTENSITY 

The  values  of  computed  from  the  many  spectra  cover  a  very  wide 

range.  We  are  then  faced  with  the  task  of  separating  and  identifying  the 
factors  responsible  for  this  variation.  If  we  consider  all  the  possible  sources 
of  influence  on  the  turbulence  intensity  at  low  levels,  we  are  led  to  the  following 

1.  Surface  roughness  (terrain  type) 

2.  Height  above  the  surface 

3.  Wind  conditions  (both  velocity  and  shear) 

4.  Atmospheric  stability 

The  Douglas  B-b6  data  provided  the  best  means  of  getting  at  the  effect 
of  different  terrain  on  low-level  turbulence.  These  measurements  were  taken 
on  specific  flight  tracks  chosen  to  represent  definite  terrain  classifications. 

The  data  taken  over  each  terrain  type  cover  a  wide  range  of  stability  and  wind 
conditions. 

The  computed  cr^.  values  from  all  the  available  spectra  were  grouped 
into  their  respective  terrain  classification  and  mean  values  calculated  for 
each  height.  A  comparison  of  the  means  for  each  group  shows  a  definite 
tendency  for  the  rms  gust  velocity  to  increase  with  increasing  surface  rough¬ 
ness.  There  is  also  some  evidence  that  the  dependence  on  terrain  rough¬ 
ness  is  a  maximum  in  the  400  to  600  ft  range  and  is  slightly  less  at  200  and 
1000  feet.  The  mean  values  of  were  then  computed,  using  all  the  avail¬ 

able  data  for  each  terrain  type.  Data  for  all  heights  within  the  terrain 
classifications  were  lumped  togethei  for  this  confutation,  since  the  small 
height  dependence  was  thought  to  be  too  indefinite  and  insignificant  to  be 
considered.  The  mean  values  falling  in  each  terrain  classification  are  given 
in  T abl  e  1 . 

The  ocean  and  water-land  flight  paths  have  a  tendency  to  have  less 
intense  turbulence  than  the  desert  classification,  even  though  the  surface 
roughness  is  not  likely  to  be  significantly  less.  The  difference  here  is 
probably  due  to  stability,  the  air  being  more  stable  over  water  where  the 


10 


AFFDL-TR-65  -  122 


water  temperature  is  generally  lower  than  the  air  temperature.  It  was 
decided  to  use  the  desert  surface  as  the  reference  and  to  reduce  all  the  #-■ 

US 

values  to  this  surface  type  by  dividing  by  the  ratio  of  the  mean  O^iox  the 
terrain  being  considered  to  the  mean  cr^  for  the  desert  terrain.  These  ratios 
are  included  in  Table  1.  The  ocean  and  water-land  types  are  considered 
to  be  in  the  desert  classification,  as  far  as  surface  roughness  is  concerned, 
a  id  hence  the  ratios  for  these  are  taken  to  be  unity.  All  the  o~u^  values  are 

reduced  by  the  appropriate  terrain  factor  Rr  from  Table  1  before  the  next 
step  in  the  analysis.  This  first  step  essentially  removes  the  effect  of  terrain 
on  the  O'  values. 


TABLE  1 

TERRAIN  FACTORS  FOR  SEVERAL  TERRAIN  CLASSIFICATIONS 


■ 

Ocean 

Water - 
Land 

De  se  rt 

harm 

V  i  rg  m 
Land 

Low 

Mins. 

High 

Mtns . 

1  .  58 

2.  32 

2.8 

3.05 

3.  20 

1 

3.  70 

4.  0 

1 

1 

1 

[ 

1 .  1 

1.15 

1 .  3 

_ i 

1 . 4 

_ 1 

According  to  the  Monin- Obukhov  similarity  theory,  the  rms  vertical 
gust  velocity  in  the  boundary  layer  for  a  given  surface  roughness  is  given  by 


(Rc  j  =  K  £  (R. ,  S) 


(10) 


where  £■  is  the  Richardson  number,  V  is  the  mean  wind  speed  at  the 
height  where  ^is  measured  and  and  6  are  universal  functions.  The 
Richardson  number  is  defined  as 


p  3  d&  ( d  V  \ 

L  e  az  \  oz  j 


(ii) 


where  O  is  the  potential  temperature.  It  can  be  shown  that 

dZ  T  \  )a  di  / 


(12) 


where  is  the  adiabatic  lapse  rate.  Combining  Equations  1)  and  12,  we 
have 
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Thus,  Equation  10  becomes 


1L) 

a  '  > di  >’ 


(13) 


From  Equation  13  we  see  that  the  quantity  o '^,/y  will  be  a  function  of 
height,  lapse  rate,  and  wind  shear  for  a  given  terrain  roughness.  The  wind 
shear  term  in  Equation  1  3  is  certainly  undesirable  since  it  is  a  parameter 
which  is  even  mere  difficult  to  obtain  than  the  lapse  rate  on  a  broad  scale. 
However,  the  wind  shear  itself  is  a  function  of  height  and  stability  along 
with  the  surface  roughness.  This  is  substantiated  by  the  tower  data  where 
a  plot  of  cr^/y  as  a  function  and  a  plot  of  /  y  versus  f?L’  show 

equally  good  correlation.  From  consideration  of  the  tower  data  it  appears 
then  that  Equation  13  may  be  rewritten  in  the  functional  form 


<^/v  ^  H(r„  +  jj,2)  (14) 

for  given  surface  roughness.  If  we  now  plot  y  as  a  function  of  height, 

assuming  the  correction  factor  RT  used  to  reduce  all  the  data  to  a  reference 
surface  is  sufficiently  accurate,  the  spread  of  data  points  at  a  given  height 
should  be  due  only  to  stability.  Figure  3  shows  the  resulting  plot  of  all  the 
available  data  including  the  tower  data.  First,  we  see  that  ^here  is  a  large 
spread  in  the  data  at  each  height.  We  also  see  that  this  spread  of  points 
increases  up  to  about  600  feet  and  then  remains  relatively  constant.  Finally 
we  see  that  the  majority  of  the  data  are  concentrated  in  a  band  of  near 
constant  cr^/R^- 

Atmospheric  stability  is  a  difficult  quantity  to  measure  and  is  certainly 
not  easily  assessed  from  flight-level  temperature  measurements  and  spot 
surface  temperatures.  The  presence  of  cumuliform  clouds  indicates  insta¬ 
bility,  but  if  the  air  is  dry,  a  superadiabatic  lapse  rate  can  develop  ir.  the 
lower  levels  without  cumulus  forming.  The  lack  of  confidence  in  the  stability 
estimates  accompanying  the  measured  spectra  is  a  major  problem  source  in 
analysis  oi  the  low-!,  vel  turbulence  spectra.  The  strong  stability  dependence 
is  evident,  but  setting  the  quantitative  stability  boundaries  without  specific 
lapse  r?te  measurements  in  the  layer  between  the  flight  level  and  the  surface 
is  difficult.  The  final  bounding  of  the  stabiiity  regions  in  Figure  3  is  based 
on  the  following  considerations; 

1.  Considering  all  the  meteorological  data  available,  it  appears  that 
the  majority  of  poix.ts  will  be  in  the  region  between  neutral  and  slightly 
unstable  lapse  rates. 

2.  The  curve  drawn  through  the  maximum  <sr/p  points  will  closely 
represent  the  line  of  maximum  possible  instability  (superadiabatic  lapse  rate) 

3.  For  maximum  stability  conditions  (inversions),  the  turbulence 
generated  in  the  boundary  layer  should  be  quickly  damped  out  with  height. 
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4.  The  tower  data  should  fit  the  over-all  picture  and  act  as  a  firm 
foundation  to  build  upon. 

5.  The  scale  of  turbulence  is  also  a  function  of  stability  and  should 
provide  an  independent  means  of  checking  the  stabilitv  boundaries . 

> 

With  the  above  items  being  carefully  considered,  the  stability  limits  are  set 
as  in  Figure  3.  The  tower  data  for  16,  40,  and  91  meters  proved  to  be  very 
helpful  in  establishing  the  lapse  rates  for  given  tr^/V  .  The  "slightly  un¬ 
stable"  line  is  drawn  through  points  for  the  tower  data  with  approximately 
dry  adiabatic  lapse  rates.  The  approximate  stability  ranges  represented 
on  F  igure  3  are:  unstable ,  3'C  per  thousand  feet ,  neutral,  1.4  Y<  3  °  C 

per  thousand  feet,  and  stable,  T  <  1 . 4°C  per  thousand  fee’.. 

The  information  indicated  by  Figure  3  relative  to  the  rms  gust  velocity 
is  summarized  below, 

1.  c increases  with  decreasing  stability  at  all  levels. 

2.  /p  decreases  with  height  under  stable  conditions. 

3.  e^./F’  remains  relatively  constant  with  height  under  neutral 
stability  conditions. 

4.  Under  unstable  conditions,  o^/v  increases  to  a  height  of  about 
600  feet  and  then  remains  nearly  constant. 

5.  At  very  low  levels,  0'ur  /  V  decreases  with  height  under  all  but  the 
most  severe  instability  conditions.  This  is  due  to  the  rapid  increase  in  V 

in  the  first  50  to  100  feet  from  the  surface  (strong  wind  shear). 


3.5  THE  SCALE  OF  TURBULENCE 

The  term  At  of  the  spectral  density  equation,  we  have  seen,  represents 
the  wavelength  of  maximum  energy.  For  the  low-level  turbulence  model 
(up  to  1000  feet)  A  is  very  nearly  equal  to  the  height.  2  under  certain  conditions. 
If  we  force  this  relationship  to  be  exact  by  letting  2 -£  for  all  conditions,  then 
the  scale  of  turbulence  becomes  A  2  where  the  scale  factor  A  ,  according  to 
the  similarity  theory,  is  a  function  only  of  stability  and  height.  All  the  '*12 
values  computed  from  the  spectra  were  divided  by  their  respective  measure¬ 
ment  heights  to  get  the  value  of  the  scale  factor  A  .  Comparison  of  the/4 
values  at  each  altitude  shows  no  consistent  relationship  to  the  associated 
terrain  type.  However,  there  does  appear  to  be  a  definite  dependence  of  A 
on  height.  The  A  values  were  plotted  against  height  only,  ignoring  the 
terrain  classification  (Figure  4).  The  existing  spread  in  data  points  at  any 
fixed  height  is  then  attributed  to  stability  variations  only.  This  result  is 
substantiated  by  the  tower  data  where  detailed  wind  and  temperature  profiles 
were  available.  The  difference  between  the  actual  lapse  rate  and  the  adiabatic 
lapse  rate,//  —  ,  was  computed  for  each  A  value.  The  correlation  between 
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the  values  and  this  liability  parameter  is  reasonably  good.  The  increase  in 
scale  parameter  A  is  approximately  linear  with  respect  to  the  stability  para¬ 
meter  (t+  The  stability  curves  drawn  on  Figure  4  were  established 

using  tne  stability  information  at  the  tower  levels,  and  extended  upward 
using  the  general  distribution  of  the  data  points  as  a  guide  in  the  same  manner 
(but  independently)  as  for  Figure  3. 

Once  these  criteria  had  been  applied  independently  to  Figures  3  and  4, 
the  stability  limits  were  adjusted  slightly  to  force  consistency  of  the  two 
figures.  Since  o-  and  A £  are  computed  from  the  same  spectra,  the  stability 
lines  should  pass  through  the  associated  points  on  the  two  diagrams.  The 
lines  are  adjusted  on  both  plots  in  such  a  way  as  to  provide  the  maximum 
correlation  between  the  data  and  yet  preserve  most  of  the  independent  stability 
analysis . 

The  curv'  s  of  Figure  4  convey  the  following  information  about  the 
behavior  of  the  scale  of  turbulence: 

1.  The  scale  of  turbulence  increases  very  rapidly  with  height  under 
unstable  conditions  in  the  boundary  layer,  and  continues  to  increase  at  a 
lower  rate  above  this  layer. 

2.  Under  neutral  conditions  the  scale  factor/?  is  nearly  constant. 

Thus,  the  increase  in  the  scale  At  is  directly  proportional  to  height. 

3.  Under  very  stable  conditions  the  scale  of  turbulence  is  nearly 
constant  with  height. 

4.  The  scale  of  turbulence  increases  with  increasing  instability  at 
all  levels . 

5.  The  scale  approaches  a  constant  value  in  neutral  and  unstable 
conditions  at  a  height  well  above  1000  feet. 


3.6  LOW-LEVEL  SPECTRUM  PREDICTION 

Figures  5  and  6  (spectrum  prediction  charts)  are  essentially  the  sanTe 
as  Figures  3  and  4  except  for  the  reduction  of  the  number  of  stability  ranges 
included.  For  practical  applications,  three  stability  ranges,  stable,  neutral, 
and  unstable,  are  considered  adequate  in  view  of  the  difficulty  in  estimating 
this  parameter.  The  solid  lines  on  these  figures  set  the  boundaries  of  the 
three  ranges  with  the  broken  lines  being  used  to  represent  each  range.  The 
stability  limits  used  for  these  ranges  are  as  follows: 

1  .  Stable  -  lapse  rate  less  than  or  equal  to  1 . 4°C  per  300  meters. 

2.  Neutral  -  lapse  rate  greater  than  1 . 4°C  per  300  meters  but  less 
than  3  °C  per  300  meters. 
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3.  Unstable  -  lapse  rate  3°C  per  300  meters  or  greater. 

Using  the  two  charts  (Figures  5  and  6)  it  is  possible  to  estimate,  with 
reasonable  accuracy,  an  expected  turbulence  spectrum  knowing  a  few  basic 
flight  parameters,  namely,  altitude,  terrain  type  (5  classifications),  expected 
stability  range,  and  the  expected  mean  wind  velocity  at  the  flight  level.  The 
scale  factor  A  is  obtained  from  Figure  6  and  multiplied  by  the  height  2  to  get 
the  turbulence  scale /4  2  .  The  factor  <ritr/R_r  (/  =  £  is  obtained  from  Figure  5  . 
Thus,  <7^  -  (e  Rr  \7)*  where  Rr  is  the  terrain  factor  and  K  is  the  predicted 
mean  wind  speed  at  Height  2  . 

3.7  GENERAL  DISTRIBUTION  OF  TURBULENCE  PARAMETERS 

The  final  phase  of  the  low-level  turbulence  model  deals  with  the 
distributions  of  the  wind  velocity  and  atmospheric  stability  over  a  broad 
area  and  subsequently  with  the  distributions  of  the  dependent  turbulence 
spectrum  parameters.  The  area  used  for  this  investigation  is  the 
continental  United  States  with  16  reporting  stations  being  used  as  a  source 
of  meteorological  data.  It  was  decided  that  the  small  differences  between 
the  farmland,  virgin  land,  and  desert-ocean  terrain  types  did  not  justify 
trying  to  classify  reporting  stations  this  carefully.  Thus,  these  three 
categories  were  combined  into  a  generally  flat  classification  with  a  terrain 
factor  of  unity.  The  exposure  ratios  for  the  three  remaining  terrain  types 
as  represented  by  the  continental  United  States  are  estimated  to  be:  flat 
terrain  55%,  low  mountains  20%,  and  high  mountains  25%.  The  reporting 
stations  were  then  chosen  to  best  represent  the  geographical  and  climato¬ 
logical  areas  of  the  United  States,  and,  at  the  same  time,  maintain  the  above 
ratios  between  the  number  of  stations  in  each  of  the  three  terrain  categories. 

The  reporting  stations  selected  are  shown  in  Table  2. 

TABLE  2 

REPORTING  STATIONS  SELECTED  FOR  METEOROLOGICAL  DATA 


Flat  (fir~f) 

Low  Mountains  { f?T  =  1. 3) 

High  Mountains  (  J?T  =  /.■?• ) 

Miami,  Fla. 
Burrwood,  La. 

Dodge  City,  Kansas 
Nantucket,  Mass, 
Green  Bay,  Wise. 
Midland,  Texas 
Oakland,  Calif. 
Bismarck,  N .  D. 
Omaha,  Neb. 

Greensboro,  N.  C. 
Seattle,  Wash. 
Pittsburgh,  Pa. 

Albuquerque,  N.  M. 

Las  Vegas,  Nev. 
Lander.  Wyo. 

Boise,  Ida. 

L 
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Two  years  of  data  covering  the  period  July  1,  1957  through  June  30, 

1959  were  processed  in  this  analysis.  All  radiosonde  and  pilot  balloon  data 
for  the  above  period  were  obtained  for  the  above  stations  on  punched  cards 
from  the  National  Weather  Record  Center  in  Asheville,  North  Carolina. 

A  computer  program  for  the  IBM  7044  was  written  to  process  this  data 
as  follows : 

1.  Compute  the  lapse  rate  in  the  layer  from  the  surface  to  300  meters 
for  each  observation. 

2.  Group  the  150  and  300  meter  wind  velocities  with  their  associateu 
lapse  rates. 

3.  Group  t.he  computed  lapse  rates  and  associated  wind  velocities  into 
the  three  stability  classes  using  the  stability  criteria  discussed  previously. 

4.  Compute  the  distribution  of  wind  velocities  at  150  and  300  meters 
associated  with  each  of  the  three  stability  classifications  for  each  of  the  four 
seasons.  Spring,  Summer,  Winte  r  and  Autumn. 

5.  Compute  the  probability  density  p(cr, for  150  and  300  meters  using 
the  wind  distributions  from  4  above  and  the  equation  cr^.  -  6  Rp/  from  the 
previous  section,  The  values  for  e  are  obtained  from  Figure  5.  The  seasonal 
grouping  is  still  maintained  here. 

6.  Combine  the  data  into  an  over-all  p  (<rw)  distribution  for  all  seasons 
for  each  of  the  stability  ranges. 

Since  only  the  150  meter  and  300  meter  winds  were  available  on  punch 
cards,  the  analysis  is  limited  to  these  two  heights  (approximately  500  and 
1000  feet) . 

The  scale  length  A Z  ,  being  independent  of  the  wind  velocity,  has 
only  a  single  value  for  each  height  and  stability  class  used. 

The  final  probability  distributions  p  are  reproduced  in  Figures  7, 

8,  and  9.  The  distributions  broken  down  into  seasons  showed  little  difference 
from  the  combined  distributions,  thus,  only  the  latter  are  given.  Figure  7 
shows  p( c 7^,)for  flight  levels  of  approximately  500  and  1000  feet  above  the 
surface  lor  stable  atmospheric  conditions.  The  distribution  for  1000  feet 
under  this  stability  category  is  concentrated  in  the  region  <r ^  less  than  0.  5 
feet  per  second.  Figure.  8  presents  the  distributions  for  heights  of  500  and 
1000  feet  under  neutral  conditions  and  Figure  9  shows  the  same  for  unstable 
conditions.  A  curious  result  is  obtained  for  the  p(crUr)  distribution  at  1000 
feet  under  unstable  conditions.  The  curve  has  a  double  peak,  one  peak 
centered  at  01^  =  4-  feet  per  second  and  the  other  at  <71^.  =  7  feet  per  second. 

The  double  peak  in  picr^)  is  a  direct  result  of  double  maxima  in  the  wind 
velocity  distribution  at.  1000  feet. 
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The  over-all  distribution  curves  shown  in  Figure  10  are  the  result 
of  combining  the  stability  categories  with  appropriate  weighting  factors.  The 
ratio  of  the  number  of  cases  in  each  stability  class  as  determined  by  the  final 
computer  output  was,  stable  43%,  neutral  36%,  and  unstable  21%.  It  was 
thought,  at  first,  that  the  radiosonde  observation  times  would  lead  to  a  bias 
tosvards  the  stable  and  neutral  categories.  However,  the  three  hour  difference 
in  local  time  across  the  country,  together  with  the  seasonal  changes  in  sun¬ 
rise  and  sunset  apparently  offset  this  bias.  An  independent  estimate  of  the 
percentage  of  time  each  of  the  stability  classes  would  apply  in  a  24-hour, 
year-round  basis,  came  very  close  to  the  percentages  above.  Therefore, 
the  composite  distributions  were  derived  using  the  computed  ratios. 

The  curves  in  Figure  10,  then,  represent  the  expected  distribution  of 
for  an  aircraft  sampling  the  whole  of  the  United  States  on  an  annual, 

24  hour  basis . 
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SECTION  4 

CONVECTIVE  STORM  TURBULENCE  MODEL 


The  main  group  of  spectra  used  in  this  analysis  came  from  the 
Douglas  B-66  turbulence  program  (Reference  9).  The  remainder  of  the 
spectra  is  taken  from  the  storm  penetration  studies  by  NASA  (Reference  10). 
There  are  88  spectra  in  all,  with  the  majority  of  penetrations  being  made 
through  the  upper  levels  in  the  range  25,  000  to  40,  000  feet.  As  in  the 
previous  low-level  analysis,  only  the  vertical  gust  spectra  are  considered. 

The  shapes  of  the  storm  spectra  are  essentially  the  same  as  those  of 
the  low-level  spectra  and  appear  to  fit  the  same  analytical  form.  There  is 
no  reason  to  suppose  that  height  has  any  relationship  to  the  scale  of  turbulence 
in  the  upper  levels  of  convective  storms,  hence  the  A  2  is  replaced  by  L  , 
and  the  storm  spectral  density  equation  becomes, 


§  (*)  = 

**  UJ 


(l-h*L)2 


(15) 


Again,  L  is  a  measure  of  the  scale  of  turbulence  and  represents  the  wave¬ 
length  of  maximum  energy  in  the  spectrum. 

The  two-point  fitting  technique  was  used  once  more  to  compute  88  pairs 
of  values  for  Cr and  L  from  the  measured  spectra.  The  results  are  plotted 
in  Figure  11.  Each  point  on  the  plot  represents  a  pair  of  data  points  ,  Li 

corresponding  to  penetration  L  .  Upon  considering  the  data  presented  on 
Figure  11,  it  is  immediately  obvious  that  a  certain  correlation  exists  between 
rms  gust  velocity  and  scale  size.  In  general,  larger  cr^'s  are  associated 
with  larger  scales.  A  closer  examination  of  Figure  11  indicates  that  there 
may  be  two  types  of  turbulence  present  in  thunderstorms;  a  low  intensity  long 
wavelength  form,  (represented  by  open  symbols)  and  a  more  intense,  shorter 
wavelength  form  (represented  by  solid  symbols).  The  two  regression  lines 
shown  in  the  diagram  were  computed  using  a  least  squares  fit  to  the  open  and 
solid  data  points  as  indicated.  The  spread  of  data  points  about  these  two  lines 
may  well  be  explained  by  the  inaccuracies  in  the  measured  spectra.  However, 
the  large  difference  in  slope  of  the  two  regression  equations  is  thought  to  be 
real.  The  unknown  parameter  which  may  be  involved  here  is  the  position  of 
the  flight  track  relative  to  the  main  updrafts  and  downdrafts.  It  is  not 
unreasonable  to  assume  that  the  most  intense,  shorter  scale,  turbulence 
exists  near  the  boundaries  of  the  main  drafts  in  areas  of  maximum  shear. 
Penetrations  made  directly  through  or  near  the  center  of  the  main  drafts  would 
lead  to  short  exposures  to  the  shear  areas  and  longer  times  spent  in  the  long 
wavelength  drafts,  possibly  accounting  for  regression  line  (1)  in  Figure  1  1. 

On  the  other  hand,  flight  tracks  along  the  edge  of  the  drafts  would  prolong  the 
exposure  to  the  highly  turbulent  boundaries  and  thus  lead  to  spectra  character¬ 
ized  by  the  regression  line  number  (2).  Unfortunately,  this  theory  cannot  be 
proven  conclusively  since  no  information  is  given  regarding  the  position  of  the 
flight  tracks  relative  to  the  active  cells  in  the  storms  penetrated. 
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Figure  II*  Distribution  of  Scale  Length  with  rms  Gust  Velocity 
For  Thunderstorm  Penetrations 
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We  can,  however,  treat  the  data  as  a  random  sampling  through  thunder- 
storms,  and  thus  consider  the  distribution  of  intensity  and  scale  parameters 
to  be  the  result  of  the  entire  sample.  The  probability  density  of  the  rms  gust 
velocity  is  computed  from  the  data  using  the  expression, 


NA  cr^ 


(16) 


where  jofar^A, 
A/  (A  <T  .  ) 

I  arc  ' 

N 


is  evaluated  at  the  center  of  the  interval  A 
is  the  number  of  points  falling  in  the  interval^ 
=  total  number  of  points 


^  =  2  feet  per  second 

The  resulting  distribution  is  shown  in  Figure  12.  The  pea  k  occurs  atao^ 
of  about  5.  5  feet  per  second  but  the  long  tail  on  the  curve  retains  a  finite 
probability  of  rms  gust  velocities  exceeding  25  feet  per  second.  This 
distribution  is  similar  to  the  distribution  given  in  the  previous  section  for 
low-level  turbulence  under  unstable  conditions  (Figure  9). 


4.  1  NUMBER  OF  GUSTS  EXCEEDING  GIVEN  MAGNITUDE 


It  has  been  demonstrated  by  Press,  Meadow's,  and  Hadlock 
(Reference  11)  that  if  turbulence  is  considered  to  be  a  stationary,  Gaussian 
random  process,  results  from  random  noise  theory,  as  derived  by  S.O.  Rice 
(Reference  12),  may  be  applied  to  gain  some  useful  information.  The 
turbulence  need  only  be  locally  Gaussian  and  stationary  over  a  short  time 
period  (a  single  penetration  through  a  convective  storm).  On  this  basis, 
an  expression  for  the  average  number  of  positive  (upward)  gusts  per  unit 
distance  may  be  w'ritten  as  follows: 

(/ 

a>  -I  '2- 


J  4  $  ur  { 4) 


No  = 


J  $us  (*) 
o 


(IV) 


Substituting  Equation  8  for  $^(4)  .  N0  can  be  integrated  in  closed  form  if 
the  upper  limit  is  finite.  In  practice,  the  aircraft  which  is  sampling  the 
turbulence  assumes  the  role  of  a  low  pass  filter  with  a  cutoff  wavelength 
w'hich  is  of  the  order  of  magnitude  of  the  wing  chord  distance.  Thus,  it 
appears  practical  to  compute  N  using  an  upper  limit  -4C  in  the  range  of  a 
few  feet  to  several  tens  of  feet.  Using  Equations  8  and  17  we  have 
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(18) 


L  >  >  1 


(19) 


For  practical  flight  altitudes  the  L  values  will  be  large  enough  to  make 
Equation  19  a  very  good  approximation.  Having  computed  the  number  of 
positive  zero  crossings,  we  can  then  go  on  to  estimate  the  expected  number 
of  upward  gusts  per  unit  distance  which  exceed  a  given  velocity  us  .  This 
may  be  expressed  as 


N  (ur)  =  M  e  3//ZO£  (20) 

o 

The  above  is  true  fora  single  value  of  Ng  and  (one  particular  turbulence 

spectrum).  It  has  been  shown  in  Reference  11  that  if  the  distribution  of  (T^j- 
is  continuous,  then  Equation  20  can  be  replaced  by 


M  (ur) 


0 


(21) 


From  the  thunderstorm  data  of  Figure  11,  we  see  that  L  has  a  distribution 
as  well  as  <7^.  Thus,  will  also  have  a  distribution.  V/e  must 

then  integrate  over  N0  as  well  as  6^.  Equation  21  then  becomes 


where  js(o^) 

■p  (No  1 <&• ) 


N 


o 


V  H  =  JfffyHZr)  dcr^  ctNo  U2) 

C  0 

probability  densi'y  of 

t  onditional  probability  density  of  Ng  for  given  cr^ 

% 

(Pc/ij  ,  number  of  positive  zero  crossings  per  unit  distance. 
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This  double  integral  was  evaluated  by  numerical  methods  using  the 
distributions  of  L  and  cr^  given  in  Figures  II  and  12.  In  the  computation 
the  actual  distribution  of  L  is  used  up  to  cr^  -  14  ft/sec.  Beyond  this  limit 
the  data  points  are  too  sparse  to  compute  a  meaningful  distribution  and  the 
number  (2)  regression  line  is  used,  is  considered  to  be  zero  beyond 

crij  =  30  ft/ sec.  The  results  of  this  computation  are  shown  in  Figure  13 
for  three  different  cutoff  wavelengths.  The  curves  are  almost  straight  lines 
of  the  form 

J-oy  M  fur)  =  M  (O)  -  rnur 

where  m  is  a  constant  and  Mfo)  depends  on  the  cutoff  waveLcngth.  Curves  for 
other  cutoff  wave  numbers  may  be  quickly  calculated  by  simply  multiplying 
the  ordinate  values  for  =  1  foot"*  by  the  square  root  of  the  value 

desired. 


NUH8ER  OF  UPWARD 


VERTICAL  GUST  VELOCITY,  us 


Figure  13.  Expected  Number  of  Upward  Gi 
Exceedina  ^  Ft/Sec  for  Thun 


AFFDL-TR-65-  122 


SECTION  5 

TURBULENCE  AT  HIGH  ALTITUDE 


In  the  portion  of  the  atmosphere  above  about  1000  feet,  we  can  distinguish 
two  classes  of  turbulence,  normally  referred  to  as  storm  and  nonstorm 
turbulence.  A  more  appropriate  designation  would  be  turbulence  in  cloud, 
referring  to  all  turbulence  in  the  troposphere  that  is  found  in  or  near  clouds, 
and  clear  air  turbulence  (CAT),  referring  to  turbulence  in  the  stratosphere 
and  in  tropospheric  regions  without  significant  cloud  formations.  Tropospheric 
turbulence  has  been  extensively  sampled  by  VGH  methods  and,  as  illustrated 
by  Section  6  of  this  report,  a  model  for  the  probability  of  occurrence  of 
specified  gust  intensities  can  be  given  with  reasonable  confidence  for  both 
storm  and  nonstorm  turbulence.  In  addition,  thunderstorm  turbulence 
spectra  are  available  from  special  observing  programs  with  instrumented 
aircraft.  A  model  for  thunderstorm  turbulence  is  given  in  Section  4, 

It  would  be  desirable  to  derive  models  for  clear  air  turbulence  and 
storm  turbulence  by  the  approach  we  have  used  in  developing  the  low  altitude 
model.  The  properties  of  the  turbulence  wouid  first  be  related  tc  meteoro¬ 
logical  factors,  and  then  climatological  data  would  be  used  to  derive  probability 
distributions  for  the  turbulence  variables.  Unfortunately,  there  are  virtually 
no  spectral  data  available  for  either  clear  air  turbulence  or  storm  turbulence 
apart  from  thunderstorms.  The  Russian  data  of  Shur  (Reference  13)  represent 
the  only  published  spectra  for  CAT.  Theory  has  not  progressed  far  enough  to 
predict  the  properties  of  the  complete  turbulence  spectrum,  or  to  definitely 
relate  turbulence  characteristics  to  meteorological  variables.  In  view  of  this 
situation,  it  i8  clearly'  impossible  to  present  any  model  for  the  turbulence 
spectrum  and  intensity  distribution  at  high  altitudes  at  this  time. 

In  the  remainder  of  this  section,  therefore,  we  restrict  our  discussion 
to  a  brief  review  of  present  knowledge  of  turbulence  structure  in  the  upper 
atmosphere,  and  discuss  the  limited  evidence  pertaining  to  spectral  shape 
and  the  meteorological  factors  that  lead  to  turbulence. 

5.  1  GENERAL  STRUCTURE  OF  THE  UPPER  ATMOSPHERE 

In  that  part  of  the  atmosphere  above  about  25,  000  feet,  observed 
turbulence  is  almost  exclusively  the  clear  air  type.  Observations  of  air  motion 
at  these  higher  altitudes  by  many  techniques  (aircraft,  balloons,  and  rockets) 
indicate  that  motions  occur  on  three  characteristic  scales.  At  the  largest 
scale  there  are  the  synoptic  circulations  depicted  on  weather  maps.  These 
have  wavelengths  in  the  range  of  one  to  several  thousand  kilometers.  The 
next  smaller  class  of  important  motions  are  mesoscale,  with  dimensions  of 
the  order  of  100  km  in  the  horizontal  and  1  km  in  the  vertical.  Between  the 
synoptic  and  meaoacales,  the  motion  of  the  upper  atmosphere  is  quite  smooth, 
ar.d  few  irregularities  exist.  The  mesoscale  motions  have  been  detected  only 
recently,  but  their  occurrence  appear#  to  be  quite  common.  Their  horizontal 
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dimensions  are  large  enough  that  they  have  limited  direct  effect  on  aircraft 
{except  when  executing  rapid  altitude  changes)  but  they  are  of  great  interest 
as  a  possible  source  of  actual  clear  air  turbulence.  Finally,  on  the  smallest 
scale,  there  is  CAT,  the  atmospheric  motions  with  scales  of  a  few  hundred 
meters  and  smaller,  CAT  occurs  less  frequently  than  the  mesoscaie 
irregularities,  though  it  is  impossible  with  current  data  to  accurately  deter¬ 
mine  where  and  how  often  it  does  occur.  Turbulence  cf  light  intensity  is 
encountered  fairly  frequently  by  aircraft,  but  it  has  been  estimated  that 
moderate  or  severe  CAT  occurs  in  less  than  one  percent  of  the  upper  troposphere 
and  lower  stratosphere  (Panofsky,  personal  communication). 

5.  2  MODELS  AND  SPECTRA  FOR  HIGH  ALTITUDE  TURBULENCE 

Numerous  indices  have  been  suggested  as  predictors  of  CAT.  Since 
these  all  represent  proposed  relationships  between  CAT  and  mean  meteor¬ 
ological  quantities,  they  can  be  thought  of  as  preliminary  models. 

Unfortunately,  none  of  the  proposed  relationships  have  been  wholly  successful, 
at  least  in  part  because  of  the  lack  of  resolution  in  meteorological  data  avail¬ 
able  for  testing  them.  Nearly  all  of  the  suggested  predictors  are  derived 
from,  or  can  be  shown  to  be  related  to,  theoretical  considerations  of  the 
stability  of  waves  or  of  the  energy  budget  of  turbulent  flow.  Wind  shear  and 
hydrostatic  stability  are  nearly  always  included.  A  number  of  criteria  for 
the  prediction  of  CAT  have  been  based  on  information  available  from  a 
weather  map  or  from  routine  aircraft  observations,  such  as  location  with 
respect  to  the  Jet  Stream,  curvature  of  the  large  scale  flow,  or  magnitude  of 
horizontal  temperature  gradients.  The  moderate  success  that  has  been 
obtained  in  relating  some  of  these  large  scale  features  to  CAT  observations 
is  almost  certainly  a  result  of  the  association  of  these  features  with  strong 
wind  shear  and  internal  frontal  zones.  It  is  suggested  that  further  improve¬ 
ment  in  models  for  relating  turbulence  to  meteorological  variables  is  unlikely 
to  .ome  from  the  introduction  of  additional  meteorological  parameters,  but 
rather  from  a  better  knowledge  of  atmospheric  fine-scale  structure  to  which 
current  theory  can  be  applied. 

At  noted  earlier,  the  only  published  spectra  for  clear  air  turbulence 
have  been  given  by  Shur  (Reference  13).  The  main  characteristics  of  these 
spectra  are  a  high  frequency  portion  that  can  be  described  by  a  -5/3  power 
law,  and  a  lower  frequency  portion  that  shows  a  steeper  slope  represented  by 
an  exponent  of  -3.  Shur  explained  the  -A  region  as  the  result  of  energy 
dissipation  by  buoyancy.  Since  Shur's  data  appeared,  Lumley  (Reference  14) 
has  presented  a  theoretical  treatment  of  the  turbulence  spectrum  for  stably 
stratified  flow  that  leads  to  the  prediction  of  a  -4  "J  range  at  low  wave  number, 
and  to  a  -A~s'3  regime  in  the  isotropic  (inertial)  sub-range.  This  substantiation 
of  Shur's  findings  is  encouraging,  but  many  more  data  will  be  required  before 
any  definite  conclusions  can  be  drawn  on  the  spectrum  of  clear  air  turbulence. 

Spectral  analyses  of  mesoscaie  motions  by  Kao  and  Woods  (Reference  15) 
have  shown  energy  decreasing  with  increasing  wave  number  more  rapidly  than 
■£~e/3.  However,  this  portion  of  the  spectrum  corresponds  to  wavelengths  in 
the  range  20-200  km,  and  is  not  reoresentative  of  the  region  considered  as 
CAT. 
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It  is  reasonable  to  expect  from  theory  and  from  the  limited  data  available 
that  CAT  spectra  should  obey  a  -k'5/3  law  for  high  wave  numbers.  (For 
engineering  application  this  region  could  probably  be  adequately  approximated 
by  a  form,  as  used  for  the  low  altitude  model.  )  Ai  lower  wave  number, 

below  perhaps  0.0015  m"'  ,  there  is  most  likely  a  steeper  region  with  slope 
exceeding  -5/3,  but  the  actual  slope  or  extent  of  this  portion  cannot  be 
determined  until  further  data  are  available. 

5.  3  POSSIBLE  PHYSICAL  MECHANISMS  FOR  THE  PRODUCTION  OF  CAT 

The  equation  for  the  rate  of  change  of  turbulent  energy  in  the  atmosphere 
shows  that  turbulence  may  be  generated  either  in  regions  of  strong  vertical 
shear  of  the  horizontal  wind,  or  in  regions  of  unstable  thermal  stratification. 
The  usual  thermal  stability  of  the  stratosphere  led  originally  to  the  notion 
that  turbulence  should  be  absent  there;  more  recently  the  frequent  observation 
of  turbulence  in  the  stratosphere  led  to  the  search  for  areas  of  unstable  or 
adiabatic  temperature  distribution.  Since  such  temperature  distributions 
have  been  shown  t  be  infrequent,  it  appears  that  CAT  is  most  frequently 
produced  by  strong  wind  shear.  Even  where  the  thermal  structure  is 
hydrostatically  stable,  strong  wind  shear  can  produce  turbulence  despite  the 
energy  required  to  overcome  the  buoyancy  forces.  The  Richardson  number, 
which  represents  the  ratio  oi  turbulent  energy  dissipated  in  working  against 
ihe  stability  to  turbulent  energy  produced  by  shear,  has  been  frequently  used 
in  attempting  to  relate  CAT  to  wind  and  temperature.  The  failure  to  achieve 
high  correlations  between  Richardson  number  and  CAT  occurrence  is  very 
probably  a  result  of  the  lack  of  ver.ical  resolution  in  upper  air  roundings. 
Computed  Richardson  numbe  rs  of  necessity  relate  to  mean  conditions  through 
an  appreciable  depth  of  atmosphere,  and  so  may  bear  little  relation  to  CAT 
produced  in  a  thin  internal  regiun  of  sharp  vertical  gradients. 

The  existence  of  very  narrow  shear  zones  which  represent  internal 
boundaries  between  different  air  masses  has  been  pointed  out  by  Danielsen 
tReference  16).  These  shear  zones  have  been  suggested  by  many  authors  to 
be  a  leading  factor  in  the  production  of  CAT.  Internal  fronts  are  characterized 
by  strong  wind  shear  and  high  stability.  The  Richardson  number  depends  on 
the  temperature  gradient  and  inversely  on  the  square  of  the  wind  shear;  thus 
the  narrower  the  frontal  layer,  the  smaller  is  the  Richardson  number  and 
the  more  likely  is  turbulence.  Typical  internal  boundaries  have  vertical 
dimensions  of  the  order  of  one  thousand  feet  or  less.  Thus  they  are  too 
narrow  to  be  depicted  in  detail  on  regular  upper  air  soundings,  but  because 
of  their  thinness,  they  are  likely  to  lead  to  small  values  of  Richardson  number. 

The  preceding  arguments  are  based  on  the  hypothesis  that  CAT  is  pro¬ 
duced  directly  through  the  growth  of  small  disturbances  in  a  favorable  large 
scale  atmospheric  environment.  An  alternative  explanation  is  that  some 
CAT  produced  indirectly  by  energy  transfers  from  mesostale  disturbances. 
These  mesoscale  circulations  have  recently  beer,  studied  by  Hildreth 
(Reference  17),  and  observations  of  such  circulations  have  been  analyzed  by 
Mantis  (Reference  18)  and  Kao  arid  Woods  (Reference  15).  It  has  been  found 
that  the  mesoscale  motions  are  definitely  not  isotropic;  both  velocities  and 
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dimensions  are  considerably  greater  in  the  horizontal  than  the  vertical 
direction.  The  exact  nature  of  the  motions  has  not  been  determined.  They 
may  be  associated  with  gravity  waves  that  result  from  perturbations  intro¬ 
duced  into  regions  of  stable  stratification  and  large  wind  shear.  Or, 
alternatively,  they  may  be  more  of  the  nature  of  turbulent  eddies,  with 
large  vorticity. 

There  is  adequate  support  in  theory  for  the  wide-spread  occurrence 
of  mesoscale  motions,  whether  they  are  waves  or  eddies.  In  Reference  17, 
it  is  shown  that  cercain  expected  combinations  of  wind  and  temperature 
structure  will  lead  to  the  growth  of  waves  from  initial  perturbations.  A 
recent  analysis  by  Pohle  et  al.  (Reference  19)  also  indicates  that  horizontal 
wind  shear  on  isentropic  surfaces  may  produce  mesoscale  motions.  It  is 
interesting  to  note  that  similar  conclusions  have  been  reached  for  tropospheric 
motions  by  Peace  (Reference  20). 

Both  the  theory  of  turbulence  (through  the  equation  for  turbulent  energy) 
arrd  classical  wave  theory  show  that  with  strong  wind  shear,  eddies  can  exist 
and  in  certain  circumstances  waves  will  be  unstable  and  will  grow.  These 
arguments  pertain  both  to  mesoscale  motions  and  to  small  scale  turbulence. 
Knowledge  of  upper  atmosphere  structure  is  adequate  for  us  to  be  certain 
that  favorable  conditions  for  both  mesoscale  and  small  scale  turbulence 
sometimes  exist;  there  remains  the  question  of  which  of  these  is  ultimately 
responsible  for  coserved  CAT. 

There  are  several  reasons  for  believing  that  mesoscale  motions, 
whether-  waves  or  vortices,  could  lead  to  CAT.  If  gravity  waves  continue  to 
grow,  they  would  eventually  "break,"  thus  leading  to  smaller  scale  turbulence. 
Eddies  on  the  mesoscale  could  also  result  in  eventual  transfer  of  energy  to 
higher  wave  numbers  and  into  small-scale  turbulence.  It  is  clear  from 
observations  however,  that  CAT  is  much  less  common  than  the  larger  scale 
motions,  so  the  problem  remains  of  determining  those  precise  conditions 
that  lead  to  CAT,  whether  it  is  generated  by  mesoscale  motions  or  directly. 

In  summary,  it  appears  that  CAT  should  be  looked  for  in  regions  where 
unstable  waves  or  mesoscale  disturbances  are  likely  to  occur,  and  in  narrow 
transition  zones  of  large  vertical  wind  shear.  Such  locations  cannot  be 
reliably  located  from  current  upper  air  data,  so  that  future  research  into 
the  causes  of  clear-  air  turbulence  should  combine  observations  or  turbulence 
with  detailed  probing  of  the  atmospheric  structure  near  the  turbulent  regions. 
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SECTION  6 

THE  EXCEEDANCE  MODEL  OF  TURBULENCE 

y 


6.  1  THE  BASIC  APPROACH 

The  presentation  of  the  atmospheric  turbulence  data  in  this  section  is 
based  on  the  theory  developed  by  S.  O.  Rice  (Reference  12)  which  enables 
one  to  calculate  the  number  of  times  that  a  random  variable,  such  as  the 
vertical  gust  velocity  ur  ,  exceeds  a  given  level.  Simple  derivations  of  this 
theory  for  a  stationary  Gaussian  (normal)  process,  the  form  which  is 
generally  considered  applicable  to  the  theory  of  the  response  of  airplanes 
in  random  atmospheric  turbulence,  are  given  in  Rice's  classic  paper  as  well 
as  in  References  21  and  22.  The  number  of  exceedances  with  positive  slope 
of  the  level  x  per  second  for  a  single  stationary  Gaussian  process  of  zero 
mean  is 


-1 


sec 


(23) 


where  x  is  the  number  of  zero  crossings  per  second  with  positive  slope, 
usually  called  the  characteristic  frequency,  and  cr  is  the  rms  value  of  x  . 
The  characteristic  frequency  may  be  calculated  from  the  spectral  density 
of  x  as  follows,  regardless  of  the  probability  density  of  x  : 
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Furthermore, 
density  of  x.  : 


the  rms  value  of  x  •  can  be  calculated  from  the  spectral 


units  of  x 


(25) 


Through  the  results  of  random  process  theory,  one  may  also  calculate  Ma  x 
and  crK  from  either  the  autocorrelation  function  corresponding  to  (to)  or 

the  joint  probability  density  function  of  x  and  x  .  In  the  present  work, 
however,  it  is  more  convenient  to  use  the  spectral  density  because  of  the 
frequency  domain  relationships  between  the  spectrum  of  the  vertical  gust 
velocity  and  the  spectra  of  the  airplane  responses. 
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The  Applicability  of  Rice's  Theory 

It  is  useful  at  this  point  to  discuss  the  applicability  of  Equations  21  and 
22  to  the  problem  of  theoretically  estimating  the  exceedance  of  an  airplane 
response  variable.  Many  studies*  for  example  References  23  and  24,  have 
shown  that  for  short  data  runs  on  the  order  of  4  minutes,  the  probability 
density  of  the  vertical  gust  velocity  is  closely  Gaussian  except  for  large 
values  of  (urza  30^,)  where  the  measured  density  is  considerably  larger 
than  predicted  by  the  Gaussian  curve  fitted  to  favor  smaller  values  of  u/~  . 
Unfortunately,  as  a  result  of  the  general  nonstationary  character  of  atmos¬ 
pheric  turbulence,  it  is  impossible  to  obtain  longer  data  runs  so  that  a  check 
of  the  probability  density  of  ur  can  be  more  reliably  obtained.  It  is  also  quite 
common  in  other  applications  that  fitted  Gaussian  density  distributions 
underestimate  the  occurrence  of  large  values  of  the  random  variable.  For 
practical  purposes,  therefore,  it  is  possible  to  regard  the  probability  density 
of  the  vertical  gust  velocity  as  being  Gaussian  for  short  periods  of  time. 

Furthermore,  it  is  significant  also  that  the  vertical  gust  velocity  has 
been  shown  to  have  an  approximately  Gaussian  character  for  a  large  range  }f 
rms  intensities  including  the  intensities  found  in  severe  thunderstorms 
(Reference  10)  (e.  g.  ,  Figure  14). 

Although  the  characteristic  frequency  N 0  x  may  be  evaluated  either 
from  measured  data  or  theoretically,  it  is  undoubtedly  the  weakest  link  in 
the  chain  of  calculations  required  for  the  theoretical  determination  of  M(x). 
This  is  understandable  in  view  of  the  large  dependence  of  N0 x  on  frequency, 
due  to  the  a)3  factor  in  the  integral  of  Equation  24.  In  fact,  because  none  of 
the  common  mathematical  forms  for  the  vertical  gust  velocity  spectrum 
approach  zero  rapidly  enough  with  increasing  high  frequency,  the  integral  tor 
the  characteristics  frequency  of  the  gust  velocity  does  not  converge.  This 
undesirable  feature  of  the  common  spectral  forms  may  be  overcome  by  making 
use  of  the  concept  of  the  "viscous  subrange"  wherein  the  spectrum  theoreti¬ 
cally  approaches  zero  amplitude  much  more  rapidly  with  frequency. 
Heisenberg's  estimate  of  this  asymptotic  behavior  is  co'7  (Reference  25). 
Theoretical  calculations  of  N0  r  for  airplane  response  variables,  although 
not  troubled  by  convergence  problems,  require  complex  linear  representations 
of  the  airplane  transfer  function,  including  consideration  of  unsteady  aero¬ 
dynamics  and  airplane  flexibility,  to  achieve  acceptable  results. 

In  a  practical  sense,  NQ  K  may  be  estimated  from  measured  data  in 
several  ways.  One  approach  suggested  in  Reference  11  makes  use  of  short 
sample  sections  of  the  time  history  of  the  variable.  Under  the  assumption 
that  the  variable  is  Gaussian  (so  that  Equation  23  applies),  one  would 
determine,  for  instance,  N  (  2  tyx  )  and  crx  from  a  short  sample  section  and 
solve  for  N0  y  using  Equation  23.  Another  approach  involves  numerical 
integrations  of  a  measured  gust  or  response  spectrum.  The  integration,  of 
course,  cannot,  be  carried  to  infinity  because  of  the  limited  frequency  range 
of  the  data,  but  the  range  of  frequencies  available  is  usually  sufficient. 
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Cumulative  Probability  Distribution  of  the  Vertical  Gust 
Velocity  for  a  Thunderstorm  Penetration  (from  Ref.  10) 
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It  may  be  assumed,  therefore,  that  Equations  23  and  24  do  in  fact  apply 
for  purposes  of  calculating  the  assumed  linear  responses  of  airplanes  to 
atmospheric  turbulence. 

The  Extension  to  Cover  Flight  Operations 

During  the  course  of  many  hours  or  miles  of  flight  operation,  an  airplane 
is  subjected  to  various  rms  levels  of  the  vertical  gust  velocity.  As  a  result  of 
this,  the  distribution  of  exceedances  or  peak  counts  is  nearly  always  much 
different  from  a  Gaussian  distribution.  In  Reference  11,  an  analytical  pro¬ 
cedure  was  developed  to  cope  with  this  situation  based  on  the  assumption  that 
the  airplane  encounters  turbulence  of  all  intensities,  that  is,  it  is  assumed 
that  the  airplane  encounters  continuous  variations  of  the  rms  vertical  gust 
velocity,  .  It  is  necessary,  therefore,  to  specify  a  probability  density 

function  and  when  this  is  done,  the  over-all  exceedance  or  peak  count 

distribution  is  given  by  the  integral 

*  !  x2  \ 

M(x)  =  A'  J  (cr^)  exp  -  —  — (26) 

*  \  U  us  *  J 


where  A „  and  Nn  v  are  independent  of  o ~t  and  o'  *=  O'2 Az  . 

X  Ot  A  1  us  X  Ur  X 

Since  p  (crr)d'<7/.represents  the  proportion  of  flight  time  that  the  airplane  spends 
in  turbulence  of  intensities  from  to  cr^  +  o( (T^  ,  the  integral  may  be 

thought  of  as  a  weighted  average  of  Gaussian  exceedance  distributions  (see 
Equation  23).  The  distribution  M  (*)  given  by  Equation  26  depends  on  flight 
condition  through  the  parameters  A '0  and  ,  which  are  functions  of  the 
scale  of  turbulence  L  and  the  airplane  transfer  function  relating  the  response 
x  to  the  gust  velocity  us  . 

Press,  Meadows  and  Hadlock  (Reference  11)  also  discuss  three  different 
probability  density  distributions  which  they  assumed  would  be  generally  useful. 
However,  in  more  recent  work  by  NASA  and  others,  only  one  of  the  three 
distributions  has  been  used:  it  is  Case  a  of  Reference  11,  which  permits 
analytical  integration  of  Equation  26.  This  probability  density  function  is  a 
Gaussian  distribution  that  is  restricted  to  positive  values  of  the  random 
variable,  and  therefore  the  amplitude  has  been  appropriately  multiplied  by  2: 

-  T-fk  ?np(-  Z  ^  50  (2,) 

Integration  of  Equation  26  with  this  probability  density  function  yields 


M[*)—  /V'  exp 
o,x 
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(28) 
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Furthermore,  it  has  been  assumed  that  there  are  at  least  two  density  functions 
(cr^j. )  to  account  for  different  types  of  turbulence.  For  example,  the  best 
known  types  of  turbulence  are  those  used  by  NASA  (e.g.  ,  References  26  and  27), 
namely,  the  nonstorm  and  storm  types  having  density  functions  and 

A  ((rur)  respectively.  The  density  functions,  both  if  the  form  of  Equation  27, 
are  combined  according  to  the  discrete  probability  of  encountering  either  one: 

=  *  Px  (Z9) 


The  probabilities  PL  ,  which  can  be  considered  the  proportions  of  total  flight 
time  spent  in  the  type  of  turbulence  described  by  -p^Cr^)  ,  are  constrained  as 
follows. 


PQ  +  P+  p  =  /  (30) 

where  P  is  the  probability  of  encountering  smooth  air  conditions.  Note  that 
since  0^**0  for  smooth  air,  pg  (&llr)*0.  Equations  exactly  analogous  to  29  and 
30  may  be  formed  if  more  than  two  distinct  types  of  turbulence  are  considered. 

For  several  distributions  with  corresponding  discrete  probabilities 

Pi  ,  Equation  28  may  be  written  in  the  form 

M(* )  ~  Na,*  £  PL  e IP  £  J  ~J^j  (31) 


which,  like  Equation  28,  is  valid  only  for  one  scale  of  turbulence  and  one  flight 
condition  (aircraft  transfer  function).  Equation  31,  when  n  is  2  or  3,  is  the 
form  used  for  the  exceedance  probability  model  in  Section  6,  except  that  the 
equation  is  put  into  a  partiy  nondimensionai  form  by  dividing  by  N 0  x  : 


The  equation  is  only  partly  nondimensionai  when  in  this  M(xj/ Na (<vs .  | X |/4X 
form  because  \x\/Ax  has  the  units  of  ft/sec. 
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Classification  of  the  Types  of  Turbulence 

For  several  reasons,  two  of  which  are  listed  below,  rather  nondesc  riptive 
names  have  been  chosen  for  the  major  types  of  turbulence.  These  are: 

P  ,  primary  turbulence, 

Pt  ,  secondary  turbulence,  and 

P,  ,  ir3  tertiary  turbulence. 

Primary  and  secondary  turbulence  may  be  associated  with  nonstorm  and  storm 
turbulence  respectively,  if  desired,  and  the  tertiary  turbulence  is  the  severe 
turbulence  associated  at  low  levels  with  certain  terrain  features  such  as 
described  in  Reference  46. 

Airplane  operations,  particularly  the  operations  of  commercial  airliners, 
have  generally  produced  peak  count  or  exceedance  data  showing  two  relatively 
distinct  distributions.  By  a  process  of  elimination  based  on  circumstantial 
evidence,  the  two  distributions  have  been  named  storm  and  nonstorm.  The 
storm  turbulence  is  identified  first  because  the  higher  intensity  distribution 
of  the  operational  peak  count  data  has  an  intensity  closely  matching  the 
intensity  measured  specially  in  thunderstorms.  The  lower  intensity  distri¬ 
bution  is  then  associated  with  nonstorm  conditions.  This  breakdown  of  the 
operational  data  is  based  on  the  assumption  that  the  rms  turbulence  velocities 
in  thunderstorms  are  distributed  according  to  Equation  27,  whereas  the  actual 
distribution  is  more  nearly  like  that  shown  in  Figure  12.  Although  the 
differences  between  Equation  27  and  the  distribution  of  Figure  12  are  large 
only  for  small  rms  intensities,  the  discrepancy  tends  to  point  out  a  need  for 
more  careful  assessment  of  the  situation. 

Another  example  of  the  arbitrariness  of  breaking  down  the  turbulence 
environment  according  to  type  of  turbulence  can  be  illustrated  by  considering 
the  fact  that  the  density  distributions  of  Cases  b  and  c  of  Reference  11 

result  in  very  good  fits  to  measured  data.  Both  of  these  distributions  yield 
exceedance  or  peak  count  distributions  exhibiting  the  same  general  character¬ 
istics  as  Equation  32.  Ordinarily,  the  measured  peak  count  and  exceedance 
distributions  curve  upward  at  large  gust  velocities  (or  accelerations)  more 
gently  than  allowed  by  the  analytical  form  of  Equation  32  with  two  or  three 
terms.  By  analyzing  the  data  available  to  them,  Press,  Meadows  and 
Hadlock  (Reference  11)  decided  that  most  peak  count  distributions  were  best 
fitted  by  the  curves  derived  using  Case  b  and  c  density  distributions  of  . 

The  data  presented  in  Figure  8  of  Reference  28  and  Figuie  6  of  Reference  29, 
Vol.  1,  foi  example,  also  show  a  more  gradual  upward  curvature.  One 
possible  interpretation  of  these  data  is  that,  instead  of  being  represented  by 
a  composite  of  2  or  3  one-sided  Gaussian  distributions  (like  Equation  27), 

0rur  is  represented  by  some  non-Gaussian  density  distribution.  It  is  evident 
that  more  research  is  necessary  in  this  area  to  better  define  the  basic 
approach. 
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The  Approach  to  the  Exceedance  Model  of  Turbulence 


The  work  described  in  this  section  deals  with  the  problem  of  determining 
the  parameters  Pf  ,  }P3  p  f of  Equation  32  as  functions  of  the 

altitude  above  the  terrain,  2  .  Tne  exceedance  probability  model  has,  there¬ 
fore,  the  functional  form 


A 1(x) 

A' 

0,K 


(33) 


The  P-  and  JlrL  parameters  are  presented  graphically  or  in  tabular  form  as 
functions  of  ?  to  define  the  basic  model.  Then  to  account  for  variation  of 
the  basic  statistics  due  to  terrain  roughness,  season,  etc-  cor rection  factors 
for  the  /7  and  ^  parameters  are  developed. 


Once  the  and  JrL  parameters  have  been  defined,  it  is  a  simple 
matter  to  obtain  the  probability  density  -fo(<Tur) and  the  cumulative  probability 
distribution  P(0^)o f  the  rms  vertical  gust  velocity  cC.  using  Equations  29 
and  27.  ” 


6.2  DETERK'  IN  AT  ION  OF  THE  P  AND  Jlr  PARAMETERS 

The  exceedance  probability  model  of  this  Section  6  was  developed 
almost  entirely  from  data  obtained  during  routine  airplane  operations  and 
during  several  special  gust  measurement  programs.  The  analysis  techniques 
varied  widely  according  to  the  nature  of  the  various  data.  Since  all  of  the 
mathematical  techniques  employed  to  analyze  these  data  require  use  of  the 
vertical  gust  velocity  spectrum,  consideration  was  given  to  the  various 
spectral  forms  that  have  been  suggested,  and  a  spectrum  form  was  selected. 
The  various  methods  used  to  determine  the  P  and  jfr-  parame:  ers  have  been 
described  in  numerous  NASA  reports,  are  not  new,  and  are  based  largely 
on  the  properties  of  Equation  32. 


Data  Sources 

Data  used  in  the  development  of  the  exceedance  probability  model  have 
been  obtained  from  many  different  airplanes  both  military  and  civil;  the  air¬ 
planes  range  in  size  from  small  fighter  airplanes  to  large  turbojet  bombers 
and  transports.  In  terms  of  the  sample  time  represented,  by  far  the  largest 
portion  is  VGH-type  data;  a  very  small  percentage  of  the  total  data  sample 
has  come  from  direct  gust  measurement  programs.  The  data  come  mostly 
from  U.S.  sources,  but  English  and  Canadian  data  are  included.  The  data 
were  collected  using  the  airplanes  listed  in  Table  3,  w'hich  is  separated 
according  to  the  type  of  data  obtained. 
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TABLE  3 

DATA  SOURCES  FOR  THE  EXCEEDANCE  MODEL 


VGH-Type  Data 


Airplane 

Reference 

Origin 

KC-  135 

30 

U.  S.  A. 

B -  52  B-F 

29 

U.S.  A. 

Sabre  Mk.  V 

31 

Canada 

DC  -4M 

32 

Canada 

Canberra  B . 6 

33 

England 

A3D 

34 

U.S.  A. 

YRB-58A 

35 

U.S.  A. 

U-2 

36,37 

U.S.  A. 

P-61 

38 

U.  S.  A. 

F  -  27 

39 

U.S.  A. 

Viscount 

40 

U.S.  A. 

Britannia 

41 

Canada 

Various  Transport 
Airplanes 

28 

U.S.  A. 

Various  Transport 

Ai rplanes 

42 

U.S.  A. 

Various  Transport 
Airplanes 

43 

England 

Direct  Gust  Data 

Airplane 

Reference 

Origin 

FH-1 

44,  45 

U.  S.  A. 

NB-66B 

9 

U.  S.  A. 

T-33 

10 

U.  S.  A. 

F-86 

10 

U.  S.  A. 

F-106A 

46 

U,  S.  A. 

B-66B 

4 

U.S.  A. 

Canberra  B, 6 

5 

England 

The  Turbulence  Spectrum  Used  for  Exceedance  Analysis 

Until  very  recently,  the  most  commonly  used  spectral  form  for  the 
transverse  components  {us,  V  )  of  isotropic  turbulence  was  the  form 
attributed  to  Dryden: 


t  (a’L> 


TC 


J+3(LSl  ) 2 
\j+(LSl )2Y 


(ft/sec)^ 

(rad/ft) 


(34) 
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which  is  here  defined  so  that 

CO 

(T^r  = J  (n,L)  ,  (ft/sec)2  (35) 

o 

This  spectrum  is  plotted  in  Figure  15  for  a  range  of  values  of  L-  .  The 
longitudinal  spectrum  corresponding  to  Equation  34  is 


(Lilf 


(36) 


Recently,  a  spectral  form  suggested  originally  by  von  Karman  has  been 
gaining  wide  approval  both  in  the  U.  S.  and  England,  The  form  exhibits  an 
__  _  S/^ 

ii  '  behavior  for  large  -f  1  and  is  therefore  consistent  with  theoretical 
results  concerning  the  change  of  spectrum  amplitude  with  frequency  in  the 
so-called  "inertial  subrange"  (References  24  and  25),  This  spectrum,  which 
is  plotted  for  various  L’s  in  Figure  16,  has  the  form 


Tur  1  ^  IT 


1  +  j  ( 1.339  l£l)2' 


137) 


for  which  Equation  35  holds.  The  cc  r  re  sponding  longitudinal  spectrum  is 


$  (nj.) 


i  ZL  _ 7 _ 

^  TC  ~[f+  '(/.3J9  Ln)jy/b 


(38) 


Certain  properties  of  the  von  Karman  spectrum,  and  the  Dryden  spectrum  as 
well,  are  given  in  Appendix  E  of  Reference  24  and  Chapter  9  of  Reference  52. 


Comparison  of  the  von  Karman  and  Dryden  Spectra 

The  most  obvious  difference  between  the  von  Karman  and  Dryden  spectra 
is,  of  course,  the  asymptotic  behavior  at  large  values  of  frequency,  the 
former  behaving  as  J2~%  and  the  latter  as  iT"2.  The  preponderance  of  recent 
experimental  evidence  favors  the  Tl"1^5  behavior,  and  this  is  in  accordance 
with  the  accepted  theory.  Although  the  von  Karman  spectrum  generally  results 
in  better  correlation  between  direct  gust  data  and  VGH-type  data  (Chapters  2 
and  10  of  Reference  52),  the  differences  between  the  von  Karman  and  Dryden 
spectra  are  small  and  considered  negligible  when  the  spectra  are  used  with 
peak  count  or  exceedance  data. 

In  addition,  the  Dryden  spectrum  facilitates  analytical  computations 
because  of  its  rational  form  (i.  e.  ,  the  numerator  and  denominator  are  simply 
polynomials  lnT-G).  Analytical  integrations  may  therefore  be  performed 
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if  desired,  and  the  spectrum  may  easily  be  factored  to  determine  a  transfer 
function  that  can  be  used  to  shape  (filter)  white  noise  to  obtain  a  time  history 
representation  of  atmospheric  turbulence  for  simulations.  Thus,  the 
rationality  of  the  Dryden  spectrum  is  a  strong  argument  in  its  favor. 

Both  the  von  Karman  and  the  Dryden  spectra  have  the  disadvantage  of 
having  too  much  curvature  near  the  knee  of  the  spectra.  The  two  spectra 
have  been  shown  to  fit  experimental  data  reasonably  well  in  the  high  frequency 
asymptotic  region,  but  the  data  shows  generally  much  less  curvature  near  the 
knee  than  do  the  spectral  forms.  In  fact,  this  is  one  of  the  reasons  that  led 
to  the  selection  of  the  simpler  spectrum  of  Equation  8  in  Section  3. 


The  Spectrum  Chosen  for  the  Present  Exceedance  Model 

In  order  to  apply  an  exceedance  probability  model  of  atmospheric 
turbulence,  a  spectrum  must  be  selected  and  corresponding  scales 
specified.  This  is  required,  first,  so  that  values  of  A/ Q  y  and  Ax  can  be 
calculated  in  order  to  determine  the  exceedance  model.  Second,  the  same 
spectrum  and  scales  used  for  obtaining  the  model  should  also  be  used  when 
applying  the  model.  However,  because  of  the  similarity  of  results  obtained 
for  exceedance  calculations  mentioned  previously,  and  because  the  turbulence 
environment  is  not  even  now  well  enough  known,  either  spectrum  may  be 
used  when  applying  the  turbulence  model. 

The  Dryden  spectrum,  Equation  34,  has  been  chosen  for  use  in  developing 
the  exceedance  model,  and  corresponding  values  of  L  are  given  in  Table  4. 
Wherever  data  were  available,  the  scales  of  turbulence  were  chosen  as  rounded 
averages  of  the  scales  determined  from  the  spectral  data  in  References  1,  4, 
5,7,  8,  9,  10,  44,  45  and  4b. 

TABLE  4 


SCALF,  OF  TURBULENCE  L  USED  FOR  DEVELOPMENT  OF 
THE  EXCEEDANCE  MODEL  OF  TURBULENCE 


Type  of  Turbulence 

~ - 

Altitude  Above  Terrain, Z,  ft 

Scale  of  Turbulence, 7_,  ft 

0  <  2  i  200 

500 

Primary 

200  <  2  <  1000 

7  50 

1000  <2<  100,000 

1000 

0  <  2  <  1000 

o 

o 

Secondary 

1000  <  2  <  10,000 

1000 

10,000  <  2  < 100,000 

2000 

Tertiary 

0  <  2  <  1000 

600 
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Reduction  and  Analysis  of  the  Data 

Equation  32  is  the  common  basis  for  all  techniques  for  evaluating 
the  P  and  J[r  parameters.  For  convenience,  the  equation  is  repeated  here 
but  without  the  term  for  tertiary  turbulence  which  will  be  treated  separately 
in  a  subsequent  section: 


M(* 

K 


=  P.  ezp 


o,x 


(-ix) 


(32) 

repeat 


The  basic  data  on  turbulence  are  almost  never  given  in  this  form,  the  usual 
forms  being  the  cumulative  frequency  distributions  of  derived  gust  velocity 
or  incremental  normal  acceleration  v,  .  These  data  are  obtained  by 
counting  normal  acceleration,  ,  peaks  from  records  and  presenting  the 
data  as  a  number  of  peak  counts  or  exceedances  per  acceleration  category 
in  a  given  altitude  band.  In  place  of,  or  in  addition  to,  the  data,  U^e 

peak  counts  may  be  given  and  the  cumulative  frequency  plots  shown.  As 
written  above,  Equation  32  may  be  considered  valid  for  the  cumulative 
frequency  distribution  of  any  airplane  response  variable  (e.  g.  )  or, 

in  fact,  for  the  gust  velocity  ur  itself.  With  each  different  variable, 
however,  the  appropriate  values  of  Noy  and  Ax  must  be  used. 

Many  reports  that  deal  with  VGH  gust  data  present  results  in  the  form 
of  cumulative  frequency  of  c/^e  per  mile,  in  the  U.  S.  A.  ,  a  standard 
procedure  has  been  used  to  convert  acceleration  peak  counts  into  Uc(e  peak 
counts.  Reference  47  describes  the  conversion  procedure  in  detaiL 
Equations  39,  40  and  41  outline  the  calculation: 


2VJ 


Acl. 


(39) 


where  the  gust  alleviation  factor 


K 

9 


is  given  in  terms  of  the  mass  rati q/^  : 


II 

0.88  A*  9 

S.ltPi 

(40) 

ZU/ 

(41) 

Usually  the  weight  iV  is  either  estimated  as  some  constant  fraction  of  the 
airplane  design  gross  weight  or  it  is  calculated  on  the  basis  of  fuel  usage. 
The  values  of  Ve  in  Equation  39  are  obtained  from  the  VGH-type  records 
for  the  same  time  at  which  the  n-y  peak  is  read. 

The  normal  acceleration  peak  count  data,  when  available,  have  been 
used  in  preference  to  U  data  ir.  the  present  work  simply  bcuaune  ul  the 
more  basic  nature  of  the  acceleration  data. 
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In  a  few  reports,  notably  References  5,  9,  10,  44,  45,  46  and  48, 
various  kinds  of  direct  vertical  gust  velocity  data  are  available.  These 
direct  data  are  considered  preferable  to  the  VGH-type  data  for  determining 
the  Jr  parameters,  but  the  P  parameters  usually  cannot  be  determined 
from  the  direct  data. 


The  Major  Source  of  Error  in  the  Reduction  of  VGH  Data 

The  acceleration  recorded  near  the  center  of  mass  of  an  airplane 
would  be  the  ideal  source  for  the  statistical  data  on  turbulence  if  the 
accelerations  were  caused  only  by  the  atmospheric  turbulence.  In  reality, 
of  course,  the  accelerations  recorded  represent  the  response  due  to 
maneuvers  as  well  as  turbulence.  The  problem  then  is  to  separate,  some¬ 
how,  the  gust  acclerations  from  the  maneuver  accelerations. 

The  approach  has  been  to  separate  the  gust  acceleration  peaks  from 
those  due  to  maneuvers  according  to  the  length  of  time  between  crossings 
of  the  reference  +1  g  level  by  the  acceleration  trace  on  the  VGH  record. 
Several  of  the  data  sources  (e.  g.  ,  References  29,  30,  35)  use  a  two-second 
interval  to  separate  maneuvers  and  gusts.  The  peak  is  called  a  gust  peak 
if  the  acceleration  trace  crosses  the  reference  level,  peaks,  and  then 
returns  to  the  reference  level  within  two  seconds.  NASA  (NACA)  has  used, 
in  the  past,  four  seconds  as  the  rule  for  the  separation  (Reference  49,  page 
■!).  However,  the  present  NASA  approach,  it  different  from  the  four-second 
rule,  is  not  described  clearly  or  with  enough  detail  in  their  reports. 

Certain  English  and  Canadian  studies  (References  31,  33,  41,  43)  have 
employed  various  counting  accelerometers  which  record,  on  counters,  the 
number  of  accelerations  greater  than  any  of  a  series  of  preselected  levels. 
Periodically  the  counters  are  photographed  or  otherwise  recorded  versus 
time.  These  counting  accelerometers  make  even  the  crudest  attempt  to 
separate  maneuvers  and  gusts  very  difficult;  in  fact,  in  some  cases,  it  is 
totally  impossible  to  effect  a  separation  of  any  kind. 

In  some  reports  in  which  it  has  been  admitted  that  separation  of 
maneuver  and  gust  accelerations  was  impossible  (Reference  50),  the  total 
acceleration  environment  has  been  presented  in  combined  form. 

Figure  17  has  been  prepared  to  compare  VGH-type  data,  as  repre¬ 
sented  by  rms-derived  equivalent  gust  velocity,  with  the  directly  measured 
gust  data,  ,  The  results  of  all  the  B-66B  low-level  gust  data  runs 
(References  48  and  51)  for  which  both  ms  values  are  given  are  plotted.  The 
rms  values  for  each  variable  were  calculated  by  integrating  measured 
spectra.  It  can  be  seen,  first,  that  the  points  are  not  centered  about  the 
line  (slope  1.  0)  of  perfect  correlation;  instead,  the  points  appear  to  be  better 
centered  about  a  line  of  slope  1.  1  that  corresponds  to  the  value  of  the  inverse 
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A 


RMS  VERTICAL  GUST  VELOCITY,  ~  FT/SEC 

Figure  17.  Comparison  of  riris  Values  of  the  Derived  Gust  Velocity,  and  the 
True  Vertical  Gust  Velocity 
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where  the  constant  C / is  the  conveision  factor  relating  to  at  a 

given  flight  condition,  with  no  pilot  input,  and  f< ^  is  defined  as 


W 


<P 


"•-5 


r  v  < 

J  -  —  (tzu)  ^  (uj)dco 

o  ar  '  ±ur 


(42) 


The  line  labeled  "trend"  in  Figure  17  represents  just  that,  a  trend  in  the  data, 
occurring  for  all  values  of  cr^  but  especially  for  the  larger  values  where 
nonlinearity  exists.  Although  they  are  not  presented  in  this  report,  similar 
data  from  the  B-66  thunderstorm  penetrations  (Reference  9)  show  the  non¬ 
linear  trend  continuing  through  substantially  larger  values  of  <Ttir  . 

Figure  17  shows,  therefore,  that  even  without  separating  gusts  from 
maneuvers,  the  normal  acceleration,  or  equivalently,  the  derived  gust 
velocity,  may  be  used  to  determine  the  parameters  for  primary  turbulence 
where  the  intensity  is  small.  This  means,  for  example,  tlv  ilues  of 

Pf  and  determined  in  Reference  50  may  be  used  to  ho1  ie  these 

parameters  for  the  exceedance  model  in  this  report. 

Finally,  results  presented  later  in  Table  6  and  Figure  ihe  data 

of  Reference  29  (two-second  rule)  and,  for  example,  for  the  data  of  References 
28  and  40  (presumably  the  four-second  rule)  show  a  large  discrepancy  which 
may  be  caused  mostly  by  the  different  crossing  times  for  separating  gusts 
from  maneuvers.  It  is  significant  that  the  average  crossing  interval  of  the 
+  1.  0  g  reference  level  for  normal  acceleration  due  to  turbulence  is  on  the 
order  of  one  second. 


Some  Examples  of  the  Determination  of  the  P  and  JP  Parameters 

Because  many  methods,  differing  only  in  detail,  were  required  to 
analyze  the  VGH  and  spectral  data  which  were  combined  into  the  exceedance 
model,  the  basic  approaches  are  described  in  a  series  of  representative 
examples  which  follow. 


The  VGH  Data  of  Reference  28 

Figure  8  of  Reference  28  (NASA  TN  D-29),  which  is  reproduced  as 
Figure  18  of  this  report,  is  a  plot  of  the  cumulative  frequency  per  mile  of 
derived  gust  velocity,  ifac  ,  peaks.  These  data  have  been  obtained  by 
combining  the  VGH  data  summarized  in  Reference  28  into  a  single  distribution 
for  each  of  several  5000-foot  intervals  of  pressure  altitude.  The  data 
represent  4.  8  x  10^  miles  of  flight  of  several  transport  airplanes  in  several 
different  airline  operations.  Because  of  the  size  of  the  combined  sample, 
great  care  was  taken  in  attempting  to  determine  the  P  and  Jr  parameters. 
However,  the  results  ot  the  analyses  made  to  determine  P  and  Jr  are 
judged  to  be  relatively  poor. 
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Specifically,  the  numerical  steps  for  determining  and  Pt  are 

described  here.  Following  the  basic  approach  given  in  Reference  27,  £  is 
calculated  according  to  the  relation 


Jr  = 
/ 


where 


C.  = 

«#,  / 


Pa 


*P 


y  is  the  inverse  slope  of  tne  cumulative  frequency  distribution. 


and 

The  calculations  of 
data: 


p^  are  based  on  the  following  assumed  average 


C-  -  4,  9t>  rad 

W  =  70,000  lb 

5  =  1200  ft2 

=  1  1.  5  ft 


The  estimated  average  true  velocity,  Kr  ,  is  given  in  Table  5.  Since  the 
calculation  of  k  ^  requires  spectral  data,  the  Dryden  spectrum  was  chosen 
with  a  scale  L  of  1000  feet  as  specified  in  a  previous  section.  The  transfer 
factor  C  was  calculated  using  the  procedure  of  Reference  47.  Moreover, 
since  the  average  data  assumed  above  do  not  warrant  detailed  calculations, 

P$  was  calculated  using  only  the  vertical  translation  degree  of  freedom  and 
was  then  multiplied  by  1.  2  to  account  for  the  pitching  and  flexibility  degrees 
of  freedom.  The  resulting  £rt  values  are  given  in  Table  5. 


The  P  values  in  Table  5  were  calculated  by  solving  the  relation 


2M(  UJe) 


(43) 


for  P1  using  the  values  of  czf  also  given  in  the  table.  Because  of  the 
close  similarity  of  the  gust  data  in  Reference  28  with  that  used  in  Reference 
27  taken  as  20  gust  peaks  (or  equivalently,  zero  crossings)  per 

mile  (page  9,  Reference  27).  The  Z’  values,  which  are  considered  poor 
at  the  higher  altitudes,  depend  greatly  on  the  values  of  <zf  which,  as  stated 
previously,  are  also  considered  poor. 

The  resulting  data,  Pf  and  ,  are  plotted  in  Figures  19  and  20 
respectively.  It  is  observed  that  at  the  lower  altitudes,  2500  and  7  500  feet, 
the  Pf  and  Jrf  values  given  in  Table  5  agree  well  with  the  final  curves  of 
Figures  19  and  20_.  At  higher  altitudes,  the  agreement  is  considerably  poorer. 
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Although  it  is  likely  that  there  is  a  large  number  of  factors  contributing  to 
these  discrepancies,  such  as  the  assumption  of  average  physical  data  for 
the  airplanes  involved  and  the  average  altitudes  and  airspeeds,  it  appears 
that  tlie  major  contributor  is  the  basic  nature  of  the  Uje  distribution  itself. 
The  distributions  shown  in  Figure  18  appear  to  contain  decreasing  amounts 
of  the  primary  turbulence  as  altitude  increases  (as  determined  by  the 
slopes).  This  is  caused  by  the  apparent  drop-off  in  the  cumulative  frequency 
of  U at  the  small  values  of  ,  which  in  t.  rn  is  caused  by  the 

distributions  of  airspeed  and  altitude  used  to  calculate  the  distributions 

from  the  y,.  peak  counts.  Hence  the  small  valups  of  Uais  .  from  which 
should  be  determined,  are  unreliable  and  clearly  cannot  be  used.  These  are 
the  points  in  Figure  18  which  are  not  connected  by  lines  to  the  corresponding 
distributions.  Curve-fitting  the  reliable  points  results  in  ^  values  which 
are  considerably  larger  than  what  is  believed  to  be  correct. 

There  are  twe  obvious  ways  to  improve  the  determination  of  £-f  and  P  . 

1.  The  simplest  is  to  work  directly  from  the  r>y  distribution, 

if  available,  arid  ignore  the  U^e  ,  and  ' 

2.  In  the  original  data  reduction  of  VGH  records,  use  a  much 

lower  reading  threshold  for  to  provide  a  reliable  estimate 

of  the  cumulative  frequency  of  Ug(a  at  smaller  values  of  U^e  . 

The  former  could  not  he  done  because  the  data  were  not  available,  and  the 
latter  is  more  in  the  form  of  a  suggestion  for  future  improvement  that 
obviously  could  not  be  effected  during  the  present  study.  A  third  way  to 
improve  the  value  of  is  described  below. 

The  VGH  Data  from  U -2  Operations  (References  36  and  37) 

The  VGH  data  for  clear  air  turbulence  at  high,  altitudes  (Reference  36) 
reported  and  supplied  (Reference  37)  by  NASA,  gives  the  total  miles  of 
flight  in  each  altitude  band  and  also  gives  the  miles  of  flight  spent  in 
turbulence.  If  the  turbulence  is  assumed  to  be  entirely  primary  turbulence, 
then  Pf  can  be  obtained  simply  by  dividing  *he  miles  of  flight  in  turbulence 
by  the  total  miles  of  flight.  Thus,  P  can  be  determined  simply  by 
inspecting  the  VGH  records. 

The  parameters  were  evaluated  as  in  the  previous  example  except 

that  was  calculated  using  a  digital  computer  program  with  an  estimated 

two-degree-of-freedom,  short  period,  representation  of  the  airplane. 

The  distributions  showed  no  recognizable  signs  of  any  secondary 

turbulence,  a  result  that  is  not  unexpected,  since  in  Reference  36  (page  4) 
it  is  stated  that,  the  operations  were  almost  entirely  in  clear  air. 
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B-66  Low-Level  Data  (References  48  and  54) 

Reference  4C  presents  turbulence  data  in  many  forms  that  would  be 
useful  for  the  development  of  the  exceedance  model,  for  example: 

1.  Peak  counts  of  the  vertical  gust  velocity,  or 

?,.  Peak  counts  of  the  derived  equivalent  gust  velocity,  U ,  and 

3,  Tabulated  spectral  variances  of  the  gust  velocity,  • 
but  only  peak  count  data  were  considered  here. 

Exceedance  analyses  of  the  B-66  results  have  been  made  by  several 
workers,  but  the  results  of  Saunders  (Reference  54)  and  Jensen  and  Hoblit 
(Reference  55)  are  more  representative  of  the  general  approach  used  here. 

They  have  achieved  similar  results  for  the  parameter  L,  using  very'  dis¬ 
similar  approaches,  Saunders  obtained  =2.72  ft/sec,  and  Jensen  and 

Hoblit:  obtained  ^  =  2.  94  ft/sec. 

Since  the  results  of  Saunders  (Reference  54),  based  on  42  representative 
data  runs,  are  considered  directly  applicable  in  the  exceedance  model 
developed  here,  a  more  detailed  discussion  is  in  order.  Saunders  obtains^ 
directly  from  the  peak  count  distribution  of  the  vertical  gust  velocity.  He 
determines  by  analyzing  the  distribution  of  the  extremes  of  peaks  from 

the  42  data  runs.  Finally,  Saunders  presents  the  constants  (  Pf  2 K'o  )  and 
(  P  2  Aj. )  (peaks/mile)  in  an  unseparated  form.  A  value  of  2^  was  there¬ 
fore  digitally  calculated  to  determine  P,  and  Pz  .  The  P  and  Jtr  parameters 
are  given  in  Table  6  for  both  primary  ani  secondary  turbulence. 


Spectral  Data  from  References  5,  44  and  45 


The  probability  density  of  cr^  given  by’  Equatior 


parameter. 

yields^ 


It  is  easily  shown  that  the  second  momer 

f 


cr 

ur 


fT  / 

jn  T  exP 


c  o nt arns  the  JC‘ ' 
the  distribution 


(44) 


The  £r  pararn  ter  may  therefore  be  considered  approximately  equal  to  the 
rms  value  of  <7^.  .  When  a  sufficiently  large  number  of  spectral  samples  is 
obtained  for  one  flight  condition,  it  is  possible  to  evaluate  b  by  simply  taking 
the  square  root  of  the  sum  of  the  spectral  variances,  CTj}.  .  The  discrete 
probability  P  cannot  be  determined  from  spectral  data. 


^This  integral  is  of  the  form  j *  2  € 

which  is  commonly  found  in  tables  of  definite  integrals.  Here  7=5^  • 

^  -j  ’  y,- 
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TABLE  6 


P  .AND  b  PARAMETERS  FOR  PRIMARY  AND  SECONDARY  TURBULENCE 


Reference 

Aircraft 

Altitude 

Z  ~Kft. 

P 

'  7 

A 

TT 

->  i 

A 

Flight 

Time'* 

hours 

30 

KC-135 

0.  5 

.  85 

2.61 

.  001 

6.  1 

47 

1.7 

.65 

2.63 

. 00066 

5.  3 

63 

3.7 

.  083 

2.  82 

39 

7.  5 

.  097 

3.  01 

-  - 

47 

15 

.  051 

2.65 

101 

25 

.  032 

2.77 

-- 

265 

35 

.  026 

2.  87 

-- 

539 

45 

.  035 

2.  32 

-- 

-- 

67 

29 

B-52 

0.5 

.  87 

3.07 

_  _ 

_  _ 

933 

0.  5 

.  67 

2.  40 

-- 

6 

1.5 

.  34 

2.63 

108 

3.  5 

.  29 

2.75 

268 

7.  5 

.  17 

2.  75 

-- 

489 

15 

.  13 

2.  41 

. 00065 

6.9 

695 

i 

25 

.  14 

2.  21 

. 00087 

5.  4 

1830 

35 

.  14 

2.  37 

-- 

7071 

42 

.  11 

2.  70 

-- 

-- 

299 

31 

Sabre  Mk.  V 

.  5 

.  703 

3.  77 

-- 

-- 

14.  1 

32 

DC-4M 

1  3 

-- 

-- 

-- 

7.  8 

1 

48.  3 

5  Canberra  B.6  .41 

-- 

3,  22 

-- 

-- 

1.04 

33  Canberra  B . 

3  .  2 

0.  40 

2.  97 

_ 

_ 

12.  34 

.  2 

0.  45 

2.  44 

-- 

-- 

5.  19 

44 

FH-1 

1 

-- 

2.  7  2 

-- 

-- 

0.  4 

45 

FH-1 

.71 

-- 

3.  83 

-- 

-- 

1.  1 
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48,  54 


B-66B 


YRB-58A  2.  5 


5.32  14.0 


36,  37 


00128 


NB-66B  27.3 


7.06  1.2 
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TABLE  6  (CONT  ’D) 


Reference 

Aire  raft 

Altitude 

Z  ~Kft. 

T"1 

A 

P 

Flight 

T  irne  'v 
hours 

40 

V  iscount 

2.  5 

.  63 

3.  70 

298 

7.  5 

.  11 

3.  54 

.  00075 

6.  43 

373 

12.  5 

.  073 

3.71 

.  00049 

7.  03 

682 

17.  5 

.  052 

2.76 

.  00051 

5.  98 

438 

22.  5 

.  030 

5.  18 

.  00030 

7. 05 

48 

28 

Transports 

2.  5 

.  153 

3.  85 

.  0035 

6.  43 

1770 

7.  5 

.  0675 

3.67 

.  002 

6.  71 

29  20 

1  2.  5 

.  0141 

4.  37 

.  0005 

8.  26 

2590 

17.  5 

.00238  4.93 

.  0003 

7.  22 

2950 

22.  5 

.00526  4.62 

.  0004 

8.  29 

1100 

43 

Transports 

1.6 

.  40  3 

2.  87 

. 00275 

6.  44 

1001 

7.  4 

.  108 

2.93 

. 00087 

5.  61 

3350 

13.  4 

.  025 

3.03 

.  000398  6. 01 

9310 

20.  1 

.  0173 

2.67 

.  000287  6. 23 

3550 

1 

33.  8 

.  0180 

2.  76 

. 00029 

3  7.  38 

2320 
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The  Parameter  P± 

The  parameter  P%  represents  the  discrete  probability  of  encountering 
the  turbulence  of  above  normal  intensity  usually  associated  with  violent  thunder¬ 
storms.  In  essence,  therefore,  determination  of  Pz  involves  the  probability 
of  encountering  a  thunder  a  turm.  What  appears  to  be  a  crude  attempt  at 
defining  Pz  was  undertaken  in  Reference  27.  At  this  time,  there  is  still  no 
method  that  wiil  assure  reliable  estimation  of  Px  as  a  function  of  height. 

All  VGH  data  obtained  from  airplane  operations,  both  civil  and  military,  are 
biased  to  an  unknown  degree  by  storm  avoidance  procedures.  Furthermore, 
there  are  no  readily  available  statistics  on  the  likelihood  of  encountering 
thunderstorms  as  a  function  of  height.  As  a  result,  the  A  curve  in  Figure 
19  has  been  determined  by  multiplying  the  Pz  values  obtained  from  the  data 
of  Table  6  by  an  appropriate  factor,  and  then  by  curve  fitting  the  modified 
P.  values.  This  factor,  approximately  equal  to  5,  was  based  on  a  discus¬ 
sion  of  the  problem  by  Notess  in  Appendix  B  of  Reference  57,  and  on  data  in 
References  58  and  59.  The  latter  two  references  show  that  storm  avoidance 
procedures  affect  essentially  only  the  value  of  P^  ,  and  not  . 

Evaluation  of  Pj  and  Aj  for  the  Low-Level  Environment 

The  tertiary  turbulence  represents  an  effort  to  incorporate  into  the 
exceedance  model  a  way  of  accounting  for  the  severe  turbulence  of  mechanical 
origin  which  is  sometimes  found  when  high  winds  blow  over  rough  terrain. 

The  estimate  Aj  -  8.7  ft/sec  is  based  straightforwardly  on  the  F-106A  High 
Intensity  Gust  Investigation  (Reference  46)  peak  count  data.  The  value  of 
Pj  -  10"^  was  estimated  by  Austin  (Reference  56)  using  data  on  the  actual 
experience  ot  B-52  fleet  airplanes  that  have  been  used  in  low-level  flight 
missions.  Austin  chose  this  value  so  that  the  F-106A  lateral  gust  velocity 

distribution  would  pass  through  the  values  of  /N0  vs.  ^'/a  that  were 

carefully  estimated  for  certain  turbulence  encounters  in  which  B-52  vertical 
stabilizers  failed.  His  basic  assumption  is  that  P  is  the  same  for  vertical 
turbulence  as  for  lateral. 

The  parameters  Pj  and  are  determined,  therefore,  for  some 

average  height  above  the  terrain  which  should  be  less  than  1000  feet.  The 
value  -  400  ft  was  selected  as  representative.  In  Figure  1  9,  fj  is 
plotted  as  a  constant  since  there  is  no  information  to  support  other  possibilities. 
Because  terrain-induced  turbulence  is  known  to  decrease  in  intensity  with 
height  above  the  terrain,  a  reasonable  but  completely  arbitrary  variation  of  £3 
with  height  is  suggested  in  Figure  20,  Note  that  no  information  is  given  for 
either  P3  or  ^  at  heights  above  1000  feet.  Use  of  the  concept  of  tertiary- 
turbulence  above  2  =  1000  feet  is  not  recommended  simply  because  of  the 
complete  lack  of  supporting  data. 


62 


AFFDL-TR-65-  122 


Compilation  of  the  Data 

The  parameters  Pf  ,  ^  ,  Px  and  lrx  are  presented  in  Table  6  for  the 
data  sources  listed  in  Table  3.  For  tertiary  turbulence  the  following  data 
were  obtained. 

P5  =  I""' 

-  8.  7  ft/ sec 

All  the  data  are  plotted  in  Figures  19  and  20  along  with  the  final  fitted  curves. 
The  curves  for  primary  and  secondary  turbulence  were  fitted  by  weighting 
the  data  points  according  to  the  following: 

1.  The  type  of  data,  VGH  or  direct  gust  velocity 

2.  The  sample  size  in  hours  of  flight 

3.  An  estimate  of  the  statistical  bias. 

The  curves  for  P  and  are  rough  estimates  and  should  be  considered 

accordingly. 

6.  3  THE  BASIC  EXCEEDANCE  MODEL  AND  ITS  USE 

The  basic  exceedance  model  is  given  by  Equation  32  which  is  repeated 
here  for  convenience: 

A  J*i\ 

A  Aj  (32) 

repeat 

in  which  the  P  and  Jp  parameters  vary  with  altitude  as  shown  in  Figures  21 
and  22,  respectively.  For  twelve  selected  altitudes  covering  the  range  from 
200  to  100,000  feet  above  the  terrain,  the  P  and  jp  parameters  are  listed 
in  Table  7  and  the  corresponding  exceedance  curves  are  plotted  in  Figure  23. 
The  exceedance  curves  are  also  called  "generalized  prediction  curves"  in 
Reference  24. 


When  the  P  and  Jir  parameters  are  specified,  it  is  possible  to  evaluate 
the  probability  density  by  using  Equation  27  along  with  an  equation  for 

the  three  types  of  turbulence  that  is  analogous  to  Equation  2  9.  Then,  further, 
an  integration  of  results  in  the  cumulative  probability  distribution 

P(<r*r)  • 

In  order  to  apply  the  exceedance  model,  for  example,  in  structural 
design  calculations,  a  spectrum  function  must  be  used.  The  Dryden  spectrum 
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TABLE  7 

EXCEEDANCE  MOjl^EL  PARAMETERS  AT  SELECTED  ALTITUDES 


Altitude 
Above 
Terrain,  Z 
f». 

Primary 

P, 

Turbulence 

A 

ft  /  s  e  c 

Secondary 

p 

Turbulence 

ft/sec 

Tertiary  Turbulence 
P  j  Jr 

-  I  ft  /  s  e  c 

200 

.72 

2.  65 

. 00016 

4. 

5 

10'5 

8.7 

1000 

.  38 

2,  90 

.  00045 

5. 

6 

10~5 

7.  0 

2000 

.  24 

3.  10 

.  0007 

6. 

n 

- 

- 

10,000 

.  07 

3.  10 

.0012 

'  • 

9 

- 

- 

20,000 

.  C7 

2.  7  5 

.  001 

8. 

9 

- 

- 

30,000 

.  085 

2.  60 

.  0007 

9. 

1 

1 

. 

40,000 

.  06 

2.  50 

.  00035 

9. 

0 

~ 

50,000 

.  02 

2.  40 

.  00011 

8. 

5 

- 

. 

60,000 

.  008 

2.  35 

2.  7  ( 1 0 )  “  5 

7. 

8 

- 

- 

70, 000 

.  002 

2.  30 

10*5 

7. 

3 

- 

- 

80,000 

.  0003 

2.  25 

2(10)"6 

6. 

4 

- 

- 

100,000 

. 00001 

2.  20 

6(10) -7 

5. 

0 

- 

- 

is  recommended  for  use  with  the  exceedance  model,  because  of  its  relative 
simplicity  when  compared  with  the  von  Karman  spectrum,  and  because  it  fits 
measured  data  nearly  as  well.  However,  the  von  Karman  spectrum  may  be 
used  if  desired;  and  furthermore,  the  Lappe  spectrum  discussed  in  Section  3 
may  also  be  used. 

The  calculation  of  bJ 0  x  depends  on  the  scale  of  turbulence  2.  ,  and 
therefore  use  of  the  values  of  L  given  in  Table  4  will  result  in  three 
separate  values  of  A '0,x  f°r  a  single  altitude.  To  realize  the  convenience 
and  simplicity  afforded  by  a  single  value  of  N Q  K  that  is  valid  for  all  types 
of  turbulence  at  one  altitude,  a  single  value  of  4-  must  be  used.  Table  8 
presents  the  recommended  single  values  of  L 
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TABLE  8 

SCALE  OF  TURBULENCE  FOR  THE  EXCEEDANCE  MODEL 


Altitude  Above  Terrain,  Z  ~  ft 

Scale  of  Turbulence,  L  — 

-  ft 

0  <  S  £  500 

500 

j  500  c.  Z  <  1000 

7  50 

j  iooo<2<  5000 

1000 

1 

j  5000c z1,  C.  100,000 

2000 

The  effect  of  neglecting  the  tertiary  turbulence  is  shown  in  Figure  23a 
for  two  pertinent  altitudes,  200  and  1000  feet  above  the  terrain.  The 
tertiary  turbulence  is  seen  to  be  important  to  the  description  of  the  environ¬ 
ment  at  heights  less  than  1000  ft.  Because  of  its  definition,  however,  the 
tertiary  turbulence  should  be  employed  only  when  airplane  operations  involve 
flight  over  rough  terrain.  Therefore,  tertiary  turbulence  would  most  likely 
be  considered  in  conjunction  with  terrain  following  operations. 


6.  4  CORRECTION  FACTORS  FOR  THE  EXCEEDANCE  MODEL 

Since  the  basic  exceedance  model  for  atmospheric  turbulence  as 
developed  to  this  point  varies  only  with  altitude  due  to  the  nature  of  the 
available  data,  an  effort  has  been  made  to  define  "correction  facto,  " ,  7n  , 
to  be  used  to  account  for  the  variation  of  the  exceedance  probabiliti  s  with 
the  following  quantities: 

1.  season  of  year, 

2.  time  oi  day, 

3.  nature  of  terrain  (surface  cover), 

4.  geographical  location,  and 

5.  terrain  roughness. 

Whereas  most  of  the  parameters  considered  for  the  spectral  model  are  of  a 
direct  meteor ological  nature,  those  listed  abi  e  are  more  of  an  operational 
nature  and,  therefore, may  be  related  to  proposed  mission  profiles  for 
airplane  s. 
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First,  assume  that  the  basic  exceedance  model  having  altitude  as  a 
parameter  represents  an  average  model,  that  is,  the  appropriate  average 
of  the  five  parameters  in  the  list.  This  model  is  represented  as 


M  (x) 


P,  ( Z)e*P 


!  Ml 


+  p 


(Z)CX-p 


!  _L  !d  ) 
I  -idK/ 


where,  as  indicated,  P  ,  Pz  ,  and  X-  are  functions  of  the  altitude 
above  the  terrain.  For  other  than  average  conditions  of  the  five  quantities, 
any  or  all  of  P  ,  ,  X-f  and  ^  may  vary.  Thus  a  given  plotted 

distribution  (e.  g.  as  in  Figure  23)  may  change  in  both  slope  and  in  its  value 
at  f X\/A  =0.  As  an  example,  over  rough  terrain  an  airplane  is  more 
likely  to  encounter  turbulence  and  the  turbulence  is  likely  to  be  more 
intense;  therefore,  both  Pt  and  should  be  larger. 


Second,  since  correction  factors  are  being  presented  here,  it  is 
logical  that  they  be  used  as  in  the  following  example. 


or 


(p) 

= 

(P) 

rough 

P 

average 

terrain 

terrain 

(i) 

- 

(Xr) 

rough 

average 

terrain 

terrain 

Furthermore,  if  corrections  are  needed  for,  say,  three  of  the  five  quantities, 
the  following  product  rule  should  be  used: 


(P) 


rough  terrain 
winter  season 
early  morning 


-  ir>  ) 

Ate: 


(mp)  _ 


time 


(  p  ! 

aver  age 
terrain, 
season, 
and  time 


Obviously,  there  are  restrictions  inherent  in  the  use  of  a  multiplication  rule 
with  the  P  va’ues.  Even  if  all  the  7r> p' s  are  greater  than  unity  simul¬ 
taneously,  P  cannot  be  greater  than  unity. 

Third,  only  gross  categories  have  been  set  up  for  evaluation  of  the 
factors.  The  data  available  at  this  time  do  not  allow  any  finely  scaled  break¬ 
down.  For  example,  season  may  be  broken  down  into  either  two  or  four 
categories.  The  two  categories  might  be  winter  and  summer,  while  the 
four  categories  might  logically  be  winter,  spring,  summer  and  fall.  Terrain 
might  be  categorized  as  in  Section  3  (5  basic  categories)  or  simply  as  Emooth 
and  rough.  The  time  of  day  may  be  broken  into  maybe  four  six-hour  periods. 
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Fourth,  some  of  the  correction  factors  vary  with  altitude  as  the 
parameters  P  ,  P}  ,  and  ^  do.  For  example,  the  effects  of  terrain 
generally  decrease  with  height  above  the  terrain. 

Finally,  along  a  different  line,  correction  factors  may  be  needed  to 
account  properly  for  the  fact  that  the  basic  exceedance  model  described 
previously  is  given  as  a  function  of  the  altitude  above  the  terrain  i  2) 
rather  than  pressure  altitude  ('  k  )  .  Unfortunately  no  data  exist  that  will 
permit  estimation  of  these  correction  factors. 


The  Correction  Factors 

Correction  factors  have  been  developed  only  for  low  altitudes  above 
the  terrain  because  it  is  only  for  these  altitudes  that  reliable  data  exist. 
Correction  factors  are  given  for  season  and  time  of  day  through  the  relation 
of  these  quantities  with  the  incident  solar  radiation  which  is  one  of  the 
primary  factors  in  convective  turbulence.  The  pertinent  data  may  be  found 
in  References  44,  51,  52,  48  and  45.  The  effects  of  season  and  time  of  day 
are  accounted  for  primarily  through  the  parameter  . 


TABLE  9 

INTENSITY  CORRECTION  FACTORS  -  TIME  OF  DAY 


j  Time  of  Day  -  Hours 

r~ 

i 

i _ 

i 

i  10  pm  to  4  am 

.  90 

1  4  am  to  10  am 

.  “  5 

10  am  to  4  pm 

1.25 

1  4  pm  to  10  r>m 

1 _ _ 

1  .  10 

TABLE  1C 

INTENSITY  CORRECTION  FACTORS  -  SEASON 


Season  (for  north 
temperate  zone) 


Months 


Oct, 

Nov, 

Dec 

1 

.  80 

,T  an, 

Feb, 

Mar 

.95 

Apr, 

May, 

J  unc 

1.  07 

J  uly , 

Aug, 

Sept 

1.  18 

7  1 
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Terrain  roughness  and  the  small  effect  of  the  nature  of  the  terrain  may  be 
accounted  for  as  done  in  Section  3  where  the  terrain  factor  R  is  defined 
from  data  in  References  48  and  51.  Here  the  symbol  f?T  is  replaced  by 
m.  but  the  data  are  the  same. 


TABLE  1 1 

INTENSITY  CORRECTION  FACTORS  -  TERRAIN  CATEGORY 


Terrain  Category 

* 

Ocean,  Water-Land,  and  Desert 

1 

F  arm 

1.  1 

Virgin  Land  (forested) 

1.  15 

Low  Mountains 

1.  3 

High  Mountains 

1.  4 

These  correction  factors  are  the  ones  for  which  consistent  trends  are  shown 
in  the  data.  However,  since  the  data  come  mostly  from  the  northern  hemi¬ 
sphere,  and  particularly  the  temperate  zones,  care  must  be  used,  and 
possibly  adjustments  made,  to  use  these  correction  factors  for  other 
locations. 

There  are  little  or  no  available  data  showing  consistent  trends  that 
might  be  used  to  determine  correction  factors  for  the  p  values. 
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SECTION  7 

CONCLUSIONS  AND  RECOMMENDATIONS 


7.  1  CONCLUSIONS 

The  Low-Level  Spectral  Model  of  Turbulence 

The  preceding  analysis  of  the  low-level  spectrum  has  shown  that  the 
rms  vertical  gust  velocity  is  a  function  of  wind  speed,  lapse  rate,  and  height 
above  the  terrain  for  a  given  surface  roughness.  This  functional  dependence 
has  been  outlined,  within  a  limited  accuracy,  for  heights  up  to  1000  feet.  The 
scale  of  turbulence  is  also  shown  to  be  mainly  a  function  of  height  and  lapse 
rate  in  the  lower  levels.  These  fairly  simple  relationships  between  the  spec¬ 
trum  and  the  meteorological  parameters  are  not  likely  to  hold  for  heights 
much  above  1000  feet. 

Similarity  theory  points  to  the  Richardson  number  as  the  important 
meteorological  factor  governing  turbulence.  At  upper  levels  beyond  the 
friction  layer,  there  is  no  longer  a  simple  relationship  between  wind  shear 
and  stability,  height  or  surface  roughness.  Thus,  any  future  work  aimed  at 
extending  this  type  of  analysis  to  higher  levels  would  very  likely  require 
knowledge  of  the  Richardson  number,  rather  than  just  lapse  rate. 

A  similar  argument  can  be  made  for  the  scale  of  turbulence.  Figure  4 
shows  that  the  scale  factor  .4  becomes  less  and  less  dependent  on  stability  and 
height  beyond  about  1500  leet.  Again,  the  parameter  which  is  likely  to  be  the 
most  useful  at  greater  heights,  where  the  wind  shear  plays  an  important  role, 
is  the  Richardson  number. 

The  Convective  Storm  Spectral  Model 

For  convective  storm  penetrations,  it  lias  been  shown  that  a  certain 
correlation  exists  between  {7  and  the  scale  length.  This  relationship  may 

be  fairly  complex  for  thunderstorm  penetration,  but  it  could  exist  in  a  simpler 
form  for  turbulence  in  clear  air  and  possibly  in  cloud  layers  other  than  cumulus 
and  cumulonimbus. 


The  Exceedance  Model  of  Turbulence 


The  exceedance  model  has  been  developed  following  the  basic  approach 
first  suggested  by  Press  and  Steiner  (Reference  27)  and  updated  by  Houbolt 
et  al  (Reference  24).  This  approach-  has  gained  wide  acceptance  because  ol 
its  simple,  yet  reasonable,  analytical  fo  m  Within  the  confines  of  this 
approach,  the  exceedance  model  developed  herein  represents  an  improvement 
over  previous  models  in  the  sense  that  a  considerably  greater  sample  of  data 
has  been  used, and  represents  an  extension,  since  the  model  is  explicitly 
defined  for  both  lower  and  higher  altitudes  than  previous  models. 
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The  model  is  presented  in  the  form  of  curves  of  the  parameters  PL 
and  Jlr^  in  Equation  32  plotted  versus  altitude  above  the  terrain.  Parameter 
P.  is  the  discrete  probability  of  encountering  the  i^h  type  of  turbulence,  and 

is  the  intensity  parameter  for  this  type  of  turbulence.  For  altitudes 
above  1000  feet,  two  types  of  turbulence,  representing  different  intensity 
distributions,  define  the  exceedance  model.  The  two  types,  called  primary 
and  seconday  turbulence,  may  be  associated  with  the  nonstorm  and  storm 
types  described  by  Press  and  Steiner  (Reference  27)  and  by  others.  For 
altitudes  of  1000  feet  and  below,  a  third  intensity  distribution,  called  tertiary 
turbulence,  is  added  to  account  for  the  severe  turbulence  generated  by  high 
winds  near  rough  terrain. 

The  parameters  P-  and  JlrL  have  been  determined  for  all  heights  above 

the  terrain  up  to  100,  000  ft.  Since  applicable  basic  data  do  not  exist  for  alti¬ 
tudes  above  approximately  70,000  ft,  extrapolations  v.ere  made.  For  all  other 
altitudes,  however,  there  were  sufficient  data  available  to  determine  the 
parameters  with  what  is  considered  adequate  statistical  reliability.  Unfor¬ 
tunately,  because  of  the  nature  of  the  data  used  in  the  exceedance  model,  it 
was  considered  impossible  to  determine  quantitative  statistical  confidence 
limits . 


The  basic  exceedance  model  varies  with  altitude  only.  Correction  factors 
have  been  developed,  where  data  were  available,  to  account  for  the  change  in 
turbulence  statistics  with  the  season,  the  time  of  day  and  with  the  nature 
(especially  roughness)  of  the  terrain.  The  correction  factors  developed 
herein  relate  only  to  changes  of  the  turbulence  intensity  parameters,  , 

since  no  consistent  data  exist  to  develop  correction  factors  for  the  P-  par¬ 
ameters  . 

The  exceedance  model  is  the  result  of  what  is  essentially  an  open-ended 
process,  that  is,  the  model  has  not  been  tested  against,  the  flight  experience  of 
an  actual  airplane  to  determine  whether  realistic  structural  design  predictions 
can  be  made.  Successful  completion  of  tests  of  this  nature  is  an  essential 
measure  of  the  reliability  of  the  exceedance  model,  as  well  as  the  spectral 
models  of  turbulence  described  in  Sections  3  and  4. 

7.2  RECOMMENDATIONS 

All  of  the  work  covered  in  Sections  3  and  4  of  this  report  has  dealt  with 
the  vertical  component  of  turbulence  only.  There  are  definite  indications, 
however,  that  the  lateral  and  possibly  the  longitudinal  spectra  may  be  des¬ 
cribed  in  a  similar  fashion.  Fortunately,  most  of  the  data  sources  used  for 
the  turbulence  models  in  Sections  3  and  4  also  contain  similar  data  for  the 
lateral  and  longitudinal  turbulence.  Primary  among  these  sources  is  the 
B-66I3  low-level  turbulence  program,  but  tower  data  and  other  airplane  data 
are  also  available.  Therefore,  since  lateral  turbulence  has  been  known  to 
cause  disastrous  failure  of  airplane  structures  at  low  levels  and  in  thunder¬ 
storms,  it  is  recommended  particularly  that  a  model  of  the  lateral  gust 
velocity  component  be  developed.  The  model  should  be  developed  so  that  it 
is  compatible  with  the  vertical  turbulence  model  in  Sections  3  and  4  of  this 
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report.  This  would  permit  use  of  the  model  of  lateral  gust  velocity  in  con¬ 
junction  with  the  vertical  gust  model. 

Separation  of  pilot- induced  normal  accelerations  from  gust-induced 
normal  accelerations  is  by  far  the  greatest  problem  connected  with  the  use 
of  VGH-type  data.  Because  the  total  current  sample  of  directly  measured 
gust  data  is  relatively  very  small,  the  VGH  data  remains  a  necessary  and 
useful  supplement  to  the  direct  data  for  developing  turbulence  models.  The 
current  practice  is  to  separate  maneuvers  from  gusts  on  the  basis  of  the  time 
interval  between  crossings  of  the  acceleration  reference  level.  A  simpler 
and  more  precise  technique  for  effecting  the  separation  can  be  obtained  by 
working  with  the  statistics,  rather  than  the  time  series  itself.  Such  a  tech¬ 
nique  requires,  however,  some  theo retical  and  experimental  research  to 
determine  for  given  flight  operations  the  pilot's  contribution  to  the  measured 
acceleration  environment.  It  is  recommended  that  this  research  be  under¬ 
taken,  not  only  to  improve  the  quality  of  gust  information  obtained  ‘rom  VGH 
data,  but  also  to  determine  how  the  interaction  between  the  gust  and  the  pilot 
affects  airplane  structural  loads  in  tu  rbulence . 

Further  improvement  of  turbulence  models  can  be  made  through  the 
universal  adherence  to  basic  standard  procedures  for  measuring  and  reducing 
turbulence  data  and  for  reporting  of  the  data.  In  view  of  current  and  future 
large-scale  programs  for  measuring  turbulence  statistics  intheU.S.A.  and 
elsewhere,  efforts  toward  standardization  are  especially  significant  and  should 
be  undertaken  as  soon  as  possible.  It  is  recommended  that  primary  emphasis 
should  be  placed  on  standards  for  direct  gust  data,  but  that  standards  for  VGH 
data  should  also  be  covered. 

Future  work  on  thunderstorm  penetrations  should  include  traverses  at 
many  different  altitudes,  m  order  to  get  some  measure  of  the  height  depen¬ 
dence  of  the  spectra.  Measurements  of  velocities  in  convective  storms 
currently  being  made  at  CAL  by  Doppler  radar  indicate  that  the  most  intense 
updrafts  are  found  in  the  upper  levels  of  the  storm  and  that,  generally,  there 
are  two  regions  of  upward  acceleration,  one  in  the  lower  levels  from  5000  to 
1  0,  000  feet  ar.d  the  other  at  mid  levels  in  the  range  ol  1  5,  000  to  20,  000  feet. 
This  would  indicate  that  the  turbulence  intensity  may  increase  upward  in  a 
similar  fashion. 

Intensity  contoured  radar  data  showing  the  flight  track  relative  to  the 
distribution  of  precipitation  in  thunderstorms  would  be  very  useful  data  to 
accompany  storm  spectra.  Analysis  of  data  available  at  present  tends  to 
show  variations  in  scale  and  intensity  which  may  well  be  related  to  the  posi¬ 
tion  of  the  f  1  i glit  track  relative  to  the  horizontal  structure  of  the  storm. 
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